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ABSTRACT 

The  motion  of  a  satellite  with  respect  to  a  fixed  system 
of  coordinates  in  space  has  been  determined.  Formulas  are  derived 
which  determine  the  aspect  of  the  satellite  axis,  and  the  aspect 
of  a  vector  perpendicular  to  the  satellite  axis.  The  telemetered 
data  consisted  of  solar  angle  measurements  in  terms  of  voltage 
from  six  sun  sensors  along  the  pitch,  yaw,  and  roll  axis,  and 
magnetic  field  measurements  from  three  mutually  perpendicular 
magnetometers.  . 

The  7094  computer  programs  to  determine  the  aspects  along 
with  the  resulting  plots  of  the  desired  angles  as  a  function  of 
flight  time  for  different  revolutions  are  exhibited.  These  plots 
showed  that  the  satellite  was  not  stable  as  expected,  but  a 
stabilizing  trend  was  noticeable  as  flight  time  increased. 
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INTRODUCTION 

To  properly  analyze  the  data  of  certain  satellite  detectors  one  must 
know  the  angle  between  the  axis  of  tics  detector  and  a  specified  vector  in 
space.  The  goal  in  this  report  is  to  determine  the  angle  between  a  detector 
perpendicular  to  the  satellite  axis  and  a  vector  from  this  detector  to  the 
earth's  center.  This  was  done  by  considering  the  angle  between  a  vector 

A 

from  the  center  of  the  earth  (which  we  will  call  U")  to  the  satellite  and 
a  vector  (which  we  will  call  e#)  perpendicular  to  the  axis  of  the  satellite 
and  in  the  opposite  direction  from  the  detector. 

In  order  to  obtain  this  angle  it  is  necessary  to  first  determine  the 
motion  of  the  satellite  about  the  pitch,  yaw,  and  roll  axes.  In  the  early 

f 

portions  of  the  flight  of  0V1-5,  motion  about  each  of  these  axes  was  active 
indicating  that  the  satellite  was  quite  unstable.  However  at  revolution 
957  the  motion  became  approximately  1.5  turns  on  the  pitch  axis  and  roll 
axis  with  rotation  of  approximately  45  degrees  on  the  yaw  axis.  At  revo¬ 
lution  1360,  complete  turns  ceased  on  each  of  the  axes  and  the  satellite 
seemed  to  be  quite  stable.  It  should  be  noted  that  a  stabilizing  trend 
was  apparent  after  revolution  957  but  the  degree  of  stabilization  could 


r 


i 


only  lie  determined  by  analysis  of  the  particular  revolution. 

The  aspect  of  the  axis  of  the  satellite  wa>  also  necessary  as  input 
data  for  the  analysis  of  the  aspect  of  the  vector  perpendicular  to  the 
axis,  that  is  the  angle  between  e^  and  U"  as  a  function  of  flight  time. 
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DESCRIPTION  OF  OV1-S  ASPECT  SYSTEM 


The  responsibility  for  designing  a  suitable  spacecraft  aspect  detection 
system  was  given  to  American  Science  and  Engineering,  Inc.  This  system 
consisted  of  six  solar  aspect  sensors  and  three  magnetometers.  For  the  • 

sun  sensors,  calibration  consists  of  determining  the  two  output  voltages 

# 

for  each  sensor  which  result  from  the  light  source  of  the  sun.  Th’e 
requirement  for  the  sun  sensors  are  a  clear  45°  conical  field  of  view  with 
axis  alignment  consistent  with  the  magnetometers.  Only  one  sun  sensor  is 
recording  at  a  given  time  and  that  sensor  is  determined  by  the  appropriate 
recorded  signature  voltage. 

A.  Position  of  Sun  Sensors  and  Magnetometers 

• 

The  0V1-5  aspect  system  sensor  locations  and  orientation  along  witlj 
the  calibration  information  is  shown  in  the  following  pages.*  The  pitch 
magnetometer  X  determines  the  component  of  the  magnetic  field  sensed  along 
the  pitch  axis;  the  yaw  magnetometer  Z  determines  the  component  of  the 
magnetic  field  sensed  along  the  yaw  axis  and  the  roll  magnetometer  Y 
determines  the  component  of  the  magnetic  field  sensed  along  the  roll  axis. 
Solar  aspect  output  #1  is  used  to  determine  the  angle  between  the  sun  and 
the  satellite  in  the  plane  of  the  reference  arrow  marked  on  the  appropriate 
sun  sensor.  Solar  aspect  output  #2  is  used  to  determine  the  angle  between 
the  sun  and  the  satellite  in  the  plane  perpendicular  to  the  arrow  marked 
on  the  appropriate  sun  sensor. 


This  information  was  provided  by  American  Science  and  Engineering,  Inc., 
Cambridge,  Massachusetts. 
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0V1-5  ASPECT  SYSTEM  SENSOR  LOCATIONS 

AND  ORIENTATION 
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4.  30 

0.  28 

0.23 

* 
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MAGNETOMETER  CALIBRATION  INI'ORMATION 


fiold  in 


Mllligauss 

X 

Y 

Z 

i  on 

4.81 

4.  80 

4.31 

550 

4.63 

4.  C2 

4.. 0  5 

500 

4.45 

4.43 

4.47 

450 

4.25 

4.  23 

4.2  5 

400 

4.04 

4.01 

4.  OS 

750 

3.  83 

3.  79 

7 .  •>  1 

300 

3.  Cl 

3.  r:7 

0  57 

250 

3.39 

3.  35 

2.34 

200 

3.  18 

3.  14 

3.13 

150 

2.98 

2.  94 

2.$3 

i 

100 

2.78 

2.  70 

1 

2.  71 

50 

2.59 

2.  53 

2.  57 

0 

2.  41 

• 

2.41 

2.4  0 

-50 

2.22 

2.  24 

2.  27 

-100 

2.03 

?.  06 

2.05 

-150 

1.84 

l.  37 

1.36 

-200 

1.63 

1.  63 

t.  07 

-250 

1.42 

1.  47 

1.  45 

-300 

1.21 

1.  26 

1.22 

-350 

.99 

• 

• 

1.03 

.  98 

-400 

* 

.78 

.  82 

• 

.75 

-150 

.  58 

.  61 

.53 

1 

C*i 

O 

o 

.38 

.  40 

3  2 

-550 

.  19 

.21 

+  .  14 

-000 

+.02 

4.  07 

-.  U3 

'S' 


I 


( 
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B.  Determination  of  Sun  Angles  from  Calibration  Curves 

Using  the  calibration  information  provided,  a  least  squares  cubic 
fit  was  made  to  each  of  the  two  output  voltages  for  all  six  sun  sensors. 
The  method  was  as  follows:  ■ 

0  =  A  +  BV  +  CV2  +  DVJ 

» 

where  V  equals  the  output  voltage,  0  equals  the  sun-sateliite  angle, 
and  A,  B,  C,  D  are  the  constants  to  be  determined.  •  * 

Let 

'n  ■  k|l  I*  *  BVk  *  CVk2  *  <  '  0kl2 


3In  3In 
3A  "  3B 


3In 

3C 


3In 

3D 


=  0 


■>An  *  Bjvk  ♦  CjVk2  ♦  Djvk3  *[0k 

♦  Bjvk2  ♦  C[Vk3  *  Djv/  .  IVkGk 
♦  CjV,4  .  DjvkS  -l\\ 


A[vk3.  Bjvk4  .  clvk5  .  Dlvk6  -  Ivk3ek 

Using  Cramer's  rule  we  can  now  find  the  desired  constants. 


n 

K 


K  K  K 


IV  Iv 


K  K 


K 


Ks 


'k 

c 

V 


Zvi5 


Zvk 


/ 


6 


lh 

ivk3 

Zvi k 

i\2  Kl 

i\4 

A  ■  Ivk2°k 

Kl  w 

i\s 

l\\ 

K  K%  ■ 

i\6 

A 

n 

l\  l\2 

i\3 

CM 

> 

- - - 

<1 

CQ 

l\\  l< 

l\s 

K\  K 

A 

n 

Ivk  K 

lvk3 

l\ 

l\2  l\\ 

V 

c  .  IVk2 

Ivk3  l\\ 

Z»k* 

• 

Zvk3 

'l\\ 

Zvk6 

.  '  A 

n 

Zvk  Zvk2 

l\ 

l\2  '  K 

lvA 

'  0  ■  i\2 

i\3  iV 

l\\ 

i\s 

z\4'  i\* 

l\\ 

A 

The  plots  of  the  different  curves  for  each  sensor  and  output  voltage  can 
be  found  in  the  appendix. 
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C.  Determination  of  the  Magnetic  Field  from  the  Calibration  Curve? 

^  i 

I ' 

A  straight  line  fit  was  made  for  each  magnetometer  using  the  appropriate  , 
calibration  information.  The  points  and  equations  used  are  shown  in  the 
appendix  along  with  the  linear  fits  to  the  calibration  curves.  It  should  be 
noted  that  the  X  magnetometer  has  its  direction  coincident  with  the  positive  * 

er  axis  in  the  previous  discussion.  The  Y  magnetometer  is  coincident  with* 

-e.,,  and  the  Z  magnetometer  is  coincident  with  -efi. 


/ 


t 
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CHAPTER  II 

DETERMINATION  OF  THE  ANGIES  BETWEEN  THE  SUN  VECTOR  AND  THE  AXES  OF  THE  SATELLITE 


A.  Theoretical  Description 


For  sensor  A  with  signature  voltage  4.69v,  we  have  the  following 


set-up: 


/ 


■A2  ♦  Z2 


/x2 


-ftl! 
.  Y2  P 


4 

V 


r  v. 

X 


Figure  2  : 

is  the  angle  determined  by  output  voltage  #1  and  02  is  the  angle  determined 
by  output  voltage  #2.  The  sun  vector  is  determined  by  the  intersection  of  the 
cones  produced  by  0^  and  a2;  and  we  will  now  determine  this  sun  vector.  However, 
it  should  be  noted  that  although  the  cones  may  intersect  in  two  places,  the 
ambiguity  is  resolved  since  we  know  which  sensor  is  reading  and  therefore  on 
which  side  the  sun  lies. 


/x2  .  z2 


cot  a 


1 


0) 


/x2  ♦  Y2 


cot  a. 


(2) 


Subtracting,  we  get 


*  2  2  2  2 
«>  X*  ♦  lc  -  cot  a 


1 

v2  v2  ,2  2 

X  ♦  Y  *2  cot  a. 


2  2  2  2  2  2 
2  -  Y  «  Y^  cot  oj  -  2  cot  a2 

t  t  t 

2“  csc“u2  *  y"  cscfcOj 
sina-, 


*Y 


sine. 


r 
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From  (1) 


v2  v2 

X  «  Y 

.  2 
cos  Oj 

2  - 
sin  a2 

-  Y2 

r  2  2i 

cos  Oj  -  sin  a2 

sin2aj 

2 

sin  aj 

.  2 

sin  a. 

9  ^  « 

' — 2 - T~ 

->X  ■*  Y  /  C0S  °1  ‘  sin  °2 


sina 


1 


Y  n  Y  Sinai 

sina. 


Z  «*Y 


sina- 

sina. 


Therefore  any  vector  R  from  the  satellite  to  the  given  light  source  can 


be  expressed  as: 


R  «  e. 


ty  /  cos^Oj  -  sin2a2l 

sina. 

1  ■ 


•e  l*y  5ln“2 


sina 


sina. 

+  %  y 


Normalizing  this  vector  and  expressing  it  in  terms  of  the  sun  vector  S 
we  get: 

-  / - 2 - 5 — 

S  ■  *ef  /  cos  Oj  -  sin  a.,  *cQ  sina2  ♦  sinaj 


1.  if  a1>0|  a2>0 


*  t  2  2 

S  ■  er  /  cos  aj  -  sin  a2  -  e0  sin|a2|  ♦  e^  sinaj 

A 

S  is  in  the  octant  ey,  -e^,  e^ 


2.  if  Oj>0j  a2<0 


S  ■  ey  /  cos  aj  -  sin  a2  -  eQ  sina2  ♦  e^  sina^ 


S  is  in  the  octant  e  .  ert,  e. 

r  o  v 


(  JV 


1 
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3„  if  c»j  <0,  a,  >0 

S  «  ey  /  cos  aj  -  sin  (kj  -  eg  sin |  a2 1  -  sin | | 

A 

S  is  the  octant  e  .  -c_,  -e. 

i  V  V 

4.  if  Oj  <0^3, <0 

S  ■  e  /  cos  Oj  -  sin  02  -  og  sina2  -  sinla^ 

* 

S  is  the  octant  e  ,  e.,  -e.  « 

r  4  4  • 

For  sun  sensor  B  with  signature  voltage  3.00v  and  output  voltages 
Bj  and  B2  corresponding  to  oij  and  Oj*  we  find 

1.  if  Bj>0,  02>o 

A  y  j  2 

S  ■-€  ♦  cos  B,  -  sin  B,  ♦  tn  sinB,  ♦  e.  sinB, 

r  l  2  0  4  4  1 

A 

S  is  in  the  ofctant  -er,  eg,  e^  | 

2.  if  Bj>0,  62<0  * 

*  /*  5  '  2 

S  »-e_  /  cos  B,  -  sin  6,  ♦  e.  sinB,  ♦  e.  sinB, 

r  1  2  0  2  4  1 

A 

S  is  in  the  octant  -e  .  -en,  e. 

r  o  ♦ 

3.  if  81<0jB2>0 

*  /  2  2 

S  *-er  /  cos  Bj  -  sin  B2  ♦  eg  sinBj  ♦  c^sinBj 

A 

S  is  in  the  octant  -e  ,  eft,  -eA 

r  o  4 

4.  ifP1<0<e2<0 

S  «-e  ,  /  cos  B,  -  sin  6,  ♦  eA  sinB,  ♦  e.  sinB, 

r  1  2  0  2  4  1 

A 

S  is  in  the  octant.  -er,  -eg,  -e^ 


/ 
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Continuing  with  this  same  approach,  we  will  now  find  S  for  sensors  C 
and  D.  for  sensor  C  with  signature  voltage  3.19v  we  have  the  following 
figure:  Y 


Z  •< 


'0 


X  r 
Figure  3 

It  should  be  noted  that  while  sensor  A  faces  along  positive  e  ,  sensor 
C  faces  along  positive  e^.  Yj  and  y?  are  the  angles  produced  by  the  output 
voltages  where  Yj  corresponds  to  output  voltige  #1. 

X2  ♦  Y2  *  Z2  cot2Y2 

Z2  ♦  Y2  -  X2  cot2Yl 


Subtracting  and  solving  for  X  we  find 

X 


*Z  sin*i 


siny, 


Substituting  (6)  into  (4)  we  find 


#22 
Y*  •  ll 


2 

cos  y- 


sin2Y, 


sin  y- 


— 5 — 

sin  y. 


(4) 

(5) 


(6) 


12 


®>  Y  - 


*  2 
siny, 


/  IT 

/  COS  Y 


2  '  Sin  Y1 


siny, 


1*1  -r 


siny, 


■>  S  =  sinyj  ♦  eQ  siny2  *  *4os^y2  -  sin^Yj 

y“2 

/cos  Yj  - 


TT 


1. 


if  Y^O,  y2>0 


S  *  *er  sinyj  +  eQ  siny 


—T~ 
sin  y2 


/ — 5 - 5~ 

S  *  er  sinyj  +  siny2  +  /cos  Yj  -  sin  y^ 


2.  if  y^O,  y2<0 


S  is  in  the  octant  e  .  ert,  e. 

roe 

A  / — g - — 7— 

S  *  er  sinyj  ♦  eQ  siny2  +  /cos  Yj  sin  y2 


3.  if  Yj<0,  Y2 > ® 


S  is  in  the  octant  e  ,  - rn ,  e. 

r  0  9 


/ — 2 - TT 

S  *  c  siny,  ♦  o„  siny,  ♦  c.  /cos  y.  -  sin 

r  1  0  2  ♦  1 


4.  if  y^O,  Yj 


<0 


S  is  in  the  octant  -e  e. 

r  0  $ 


/ — 2 - 

S  ■  e^  sinyj  ♦  eQ  siny2  ♦  /cos  y^  -  sin  y2 


S  is  in  the  octant  -e  ,  -e_,  e. 

r  0  9 


For  sensor  D  with  output  signature  voltage  2.42  and  output  angles 
61  and  62  respectively,  the  sun  vector  can  now.be  expressed  as 


S  *  e  sinfi.  -e.  sin6,  -  e. 

.  T  10  2  $ 


/ — 3 - 5 — 

/cos  6j  -  sin  62 


This  follows  from  the  analysis  for  sensor  C  except  that  eQ  is  replaced 
by  -eQ  and  e^  by  -e^. 

The  reader  can  determine  which  quadrants  the  sun  vector  lies  for  (8) 
depending  upon  the  conditions  placed  on  6^  and  62> 


*  (7) 


(8) 
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Now,  wo  will  find  S  for  sensors  E  and  F. 
voltage  ,86v,  we  have  the  following  figure: 

Y 


Sensor  E  faces  positive  e0  and  5^  and  £2  are 

the  respective  output  voltages. 

7  2  2  2 

*  l*  -  r  cot*  Cj 

2  2  2  2 

+  z *  r  COt^  52 

Subtracting  and  solving  for  Y  we  get 

Y  =*  X  Sing? 
sintj 

Substituting  (10)  into  (9)  ,  we  find 

z2  =  x2  (cos2 5 j  -  sin2C2) 

- 2 - 

sin  Cj 

z  *  X  *4os2C.  -  sin2C- 
sin^i 

X  =  X  sinCj 

sin^ 


For  sensor  E  with  signature 


These  cones  actually  intersect 
as  can  be  seen  from  the  figure 

for  sensor  C  if  we  let  e\  =-e,., 

<t>  .  O’ 

e  =  e.,  e  =  e  where  the  bars 
0  <r  r  r 

represent  the  unit  vgctcfrs  for 
sensor  E  (and  then  remove  the 
bars) . 

the  angles  corresponding  to 

(9) 

(10) 


/ 
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The  normalized  sun  vector  can  now  be  written  as 

,  /  2  2 

S  =  er  sin^j  *  cQ  /  cos  ^  -  sin  £2  1  e$  sin£2 

Due  to  output  voltage  #2  of  sensor  E  as  can  be  seen  by  the 

calibration  of  output  voltage  #2, 

-  1 - > - T~ 

S  *  e^  sm£.  +  e_  /  cos  £.  -  sin  £„.  +  e.  sin?- 

Once  again  the  reader  can  determine  which  quadrants  the  sun  vector  lies 
in  for  (11)  depending  upon  the  conditions  placed  on  £j  and  £2* 

For  sensor  F  with  output  signature  voltage  1.57v  and  output  voltages 
fj  and  f2  respectively. 

a  /  2  2 

S  =  e  sxnf,  -  en  /cos  f,  -  sin  f,  +  e.  sinf_ 

r  1  0  1  2  v  2 

This  follows  from  the  analysis  for  sensor  E  except  that 

eQ  is  replaced  by  -eQ  and  e^  is  not  changed  due  to  the  remark 
preceding  (11)  . 

A 

In  summary,  for  each  of  the  sun  sensors,  the  unit  vector  S  can  be 
expressed  as  follows: 

A 

Sensor  A  S 

A 

Sensor  B  S 

A 

Sensor  C  S 


Sensor  D  S 

A 

Sensor  E  S 


Sensor  F  S 


/  2  2 

*  e  /  cos  a,  -  sin  a_  -  e,,  sina_  +  e.  sina. 

r  1  2  0  2  v  1 


-e  /— 2-  — 2 


cos  6,  -  sin  ♦  e.  sinB,  +  e.  sinB. 
r  1  2  0  2  9  1 


*  er  sinYj  ♦  e0  siny2  +  /  cos2Yl  - 


.  2 

sin  y , 


/2  2 
cos  6.  -  sin  6, 
r  i  u  i  9  l  l 


/2  2 

r  ^  cos  £j  -  sin  C.2  +  e^  sin£. 


/2  2 
cos  f,  -  sin  f„  +  e.  sinf» 


kl  0 


2  *  2 


(ID 


(12) 


/ 
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To  determine  the  angle  between  the  sun  vector  and  the  axes  of  the 
satellite,  we  need  only  consider  the  respective  dot  products.  That  is: 

S.er  =  cos  0s,  ef)  =  cosine  of  the  angle  between 

S  and  e  . 
r 

S.eQ  *  cos  0s,  e0)  =  cosine  of  the  angle  between 

S  and  Cq. 

A 

S.e^  =  cos ()S,  e^)  =  cosine  of  the  angle  between 

S  and  c . . 

9 


* 


Using  this  approadt,  we  find  the  following  angles: 


Sensor  A 

Sensor  B 

Sensor  C 

A 

OS,  e,r)  arcos  * 

/  2  2 
'  cos  Oj  -  sin 

/  2  2 
arcos  (-  /  cos  ctj  -  sin  a2  ) 

90-Yj 

| 

OS,  e0) 

90  +  a2 

90- e2 

90-y2  I 

(*S.  e4) 

90  -•  aj 

90-Bj 

arcos  / 
• 

'  2  2 
cos  Yj  •  sin  Y2 

Sensor  D 

Sensor  E 

Sensor  F 

OS,  er) 

90-fij 

90*62 

arcos (- 

/  i  2 

/  cos  6j  -  sin  6 

OS,  e0) 

90-Cj 

/  5  T 

arcos  /  cos  -  sin  52 

90-C2 

OS,  e#) 

90- fj 

. .  . . — 

arcos (-  /  cos  fj  -  sin  f2  ) 

90- f 2 

In  each  case  output  voltage  #1  corresponds  to  the  angle  subscripted  with 
a"l"  and  similarly  for  output  voltage  #2. 

4 


•« 
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B.  Results  and  Plots  for  Different  Revolutions  (Full  Orbits  and  Real  Time) 

As  can  be  seen  by  the  following  plots,  the  satellite  is  not  well-behaved. 
During  the  later  orbits,  the  motion  of  the  satellite  seems  to  becoming  more 
stable  than  that  of  Rev.  480,  however,  at  no  time  can  a  definite  precession 
angle  be  found.  The  X  axis  of  each  of  the  plots  represents  seconds  Greenwich 
Meridian  time.  and  the.  Y  axis  represents  the  angle  in  degrees.  * 

The  program  to  determine  the  angles  between  the  sun  and  the  axes  of  thjj 

i 

satellite  is  incorporated  into  the  major  program  that  will  be  discussed  in 
Chapter  V.  The  angles  that  are  found  range  from  0°  to  180°,  so  to  predict 
a  complete  turn  on  either  the  pitch,  roll,  or  yaw  axis  one  would  have  to 
examine  the  appropriate  plot.  An  example  of  this  appears  in  the  plot  of 
()tS,  e^)  for  revolution  480.  The  sharp  descent  at  approximately  53.4k  seconds 
indicates  a  complete  turn  of  the  roll  axis,  i.e.  ()S,  e^)  is  going  from 
positive  0"  -*■  180*  to  negative  180*  -*  0*. 

Since  the  signature  voltages  for  the  sun  sensors  were  not  extremely 
accurate,  limits  were  set  on  each  signature  voltage  to  determine  the  appropri¬ 
ate  sun  sensor  reading.  These  limits  and  the  action  taken  can  be  found  in 
the  program  referenced  above. 
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CHAPTER  III 

DESCRIPTION  OF  THE  FIXED  REFERENCE  SYSTEM 
A.  Fixed  Reference  System  with  Respect  to  the  Vernal  Equinox 

Since  we  are  considering  the  motion  of  a  satellite,  the  contribution 
of  the  earth's  rotation  about  its  axis  and  its  rotation  about  the  sun  must 
be  taken  into  account  when  trying  to  determine  the  aspect  of  the  satellite, 
with  respect  to  a  fixed  system  of  coordinates.  To  express  the  aspect  in 
terms  of  a  gaocentric  system  of  coordinates  would  not  have  much  meaning4 
due  to  the  above  angular  contributions.  The  fixed  system  used  is  with 
respect  to  the  vemal  equinox*  (March  21) . 

Let  i  be  a  unit  vector  parallel  to  the  line  from  the  observer  to  the 
point  on  the  celestial  sphere  where  the  sun  appears  at  the  time  of  the 
vernal  equinox  and  directed  toward  £he  sun.  Let  k  be  a  unit  vector  parallel 
to  the  polar  axis  of  the  earth  and  j  *  kxi 


Figure  17 

The  unit  vector  j  lies  in  the  equatorial  ,-lane  of  the  earth. 

Every  vector  in  this  fixed  system  can  bo  expressed  in  the  form 

A 

V  C0SQV  C0S*v  +  j  cosOy  sin*y  *  k  sin0y  . 

V  is  an  arbitrary  unit  vector  in  this  fixed  system,  %  is  the  angle  between  the 

*  *  4 

equatorial  plane  and  the  vector  V,  is  the  azimuth  of  V  with  respect  to  the 
vemal  equinox. 

*The  vemal  equinox  is  defined  as  the  intersection  of  the  equatorial  plane  of 
the  earth  and  orbital  plane  of  the  earth. 


/ 

! 
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B.  Expression  of  Required  Vectors  in  this  Fixed  System  of  Coordinates 


In  this  section  we  will  describe  the  unit  vectors  S,  M,  e  ,  U",  and  e.  in 

T  *  9 


the  fixed  system  of  base  vectors  i,  j,  k  that  will  be  uised  in  the  next  three  chapters.1 

1.  The  sun  vector 


The  unit  sun  vector  S  from  the  earth  to  the  sun  can  be  expressed 


as 


S=  i  cosO,  cos*  ♦  j  cosO  sin*  +  k  sino. 

s  s  s  s  s 


(13) 


where  0  is  the  apparent  declination  of  the  sun  and  ♦  is  the  apparent 

S  Si 


right  ascension  with  respect  to  the  vernal  equinox  at  March  21.  These 
angles  were  found  in  The  American  Ephemeris  and  Nautical  Almanac  -  1966. 

2.  The  magnetic  field  vector  , 


The  unit  vector  M  representing  the  magnetic  field  will  be  expressed 


as 


i  cos0H  cos*H  ♦  j  cos0H  sin^  +  1c  sin0H 


The  determination  of  these  angles  0^  and  *H  for  the  fixed  system , of 


coordinates  i,  j,  k  will  be  discussed  in  chapter  IV. 

3.  The  axis  of  the  satellite,  e_  . 


The  axis  of  the  sateEite  can  be  expressed  as  the  unit  vector 


e_  * 


i  cos0  cos*  +  j  cos0  sin*  +  k  sine 

”  i 

In  Chapter  V,  we  will  go  into  a  detailed  description  explaining 


how  to  obtain  the  angles  0  and  *. 

A 

4.  The  vector  U"  from  the  center  of  the  earth  to  the  satellite 


Given  the  following  information: 
tm  »  eppemeris  transit  time 

|  • 

0  ■  latitude  of  the  satellite  from  the  ephemeris 
*  ■  longitude  qf  the  satellite  from  the  ephemeris 
t  ■  Greenwich  Meridian  time  in  secs. 


*s  ■  right  ascension 


0  ■  angle  between  the  nose  axis  of  the  satellite  and  the  equatorial 


plane  ' 

*  ■  azimuth  of  the  nose  axis  W.R.t.  the  vernal  equinox 


(14) 


(IS), 


i  i 


.  i 


t 

i, 

Y: 


f  ■ 


I 


I 


I 


I 
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2tl  '  /  « 

w  «  vT  where  T  is  secs,  in  a  sidereal  day 

/  '  1  ) 

.  i  1  * 

Given  the  angles  0£  and  we  can  set  the  bnit  vector  U"  in  a 
rotating  system  as: 

i  A  .  •• 

U"  =  i  cos0£|  cps4>E  +  j  eosO  sin4»E  ♦  k  sinO£ 

where  i  is  in  the  Greenwich  Meridian  Plane,  k  is  in  the  direction  of  the 

north  pole  and  j  ■  k  x  i.  Both  i  and  j  lie  in  the  equatorial  plane  of 

the  earth.  .  1  i 

If  wc  let  tm  equal  the  ephemeris  transit  tims  at  which  the  Creenwich  « 

1  Meridibn  Plane  transits  the  sun  line  (this  time  ican  be  found  on  pages  19-33 

in  the  above!  mentioned  almanac),  and  us  =  where  T  is  the  time  in  seconds 

i  * 

for  a  sidereal  day,  and  t  >  GMT  in  seconds,  then 
'  '  '  t  i  '  ' 

I  ■  i  cos J(D(t-tm)  +  *s]  +  j  sinfu(t-tm)  +  $s] 

i  1 

3  »-i  sin[(<)(t-tm)  ♦  t  ]  +  j  cos[w(t-tm)  +  *„] 

I  t  S  S' 

Using  1  (17)  and  (18)  we  can  now  express  U"  in  the  i,  j,  k  fixed  ! 

•  j  i  ) 

system,  i.e. 

'  ' 

U"  ■  i  cosCjj  cos[u(t-tm)  ♦  *s  ♦  *g]  +  3  cosOg  sin[w(t-tm)  +  4>s  ♦  $E] 

i  1  { 

i  *  k  sin0B 

I 

where'  k  is  parallel  to  it.  Perhaps  at  this  point,  some  diagrams  might  help 


(16) 


(17) 

(18) 


(19) 


'  Figtire  16 


;  > 


*  A  sidereal  day  is  the  duration  of  one  rotation  of  the  earth  on  its  axis  with 
respect  to  the  vernal  equinox.  A  sidereal  day  is  23  hours,  56  minutes,  4.09054  secs, 
of  mean  solar  time.  I  < 


!  I 


■  I 


/ 
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In  this  diagraa  j  ■  k  x  i, 
along  with  i.  Looking  at  tho 


and  j  lies  in  the  equatorial  plane  of  the  earth 
above  figure  froa  another  view  wc  have 


5.  The  unit  vector  et  in  the  fixed- system. 

Given  the  anglis  6  and  t  any  vector  R  along  the  nose  axis  of  the  satellite 
in  the  fixed  systea  can  be  written  as 


R  •  r(i  cosO  cost  ♦  j  cost)  tint  ♦  k  sinO) 

3R 

■>er  ■  ■  i  cosG  cost  ♦  j  cos©  sint  ♦  k  sin0 

*  1  an 

Nj  »  -  *§0  sinO  cos#  -  j  sinO  kin#  ♦  k  cos@ 

n2  “  rb?  U  -i  sin#  * J  co#* 

and  «r,  hj,  N2  defines  an  orthogonal  system. 
Since  the  satellite  aay  be  spinning,  e^  is  in  the  plane  of 

a  mm 

and  and  e^  ■  Nj  cosp  ♦  N2  sinp. 


(20) 

(21) 


>1 


I*S  #» 


er  nose  axis 


Figure  20 


It  should  be  noted  that  we  could  have  let  p'  ■  angle  between  N 2  and  e#. 
Then  we  would  get 


e#  ■  Nj  sinp'  ♦  N2  cosp* 


Using  the  angle  p  we  get 


e#  «  Nj  cosp  4  N2  sinp  ■  (-i  sine  cos*  -  j  sine  sin*  4  k  cose)  cosp 

4  sinp(-i  sin*  4  j  cos$) 

■>e#  «-i(sin0  cos*  cosp  4  sin*  sinp)  -  j (sine  sin*  cosp  -  cos*  sinp) 


4  k  cose  cosp 


All  that  is  needed  to  uniquely  determine  this  unit  vector  is  the  angle 
p,  and  this  will  be  discussed  in  Chapter  VI. 
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CHAPTER  IV 

DETERMINATION  OF  AND  $H  FOR  A  FIXED  SYSTEM 

In  this  chapter  we  will  determine  the  magnetic  field  of  the  earth  in 
a  fixed  system  of  coordinates.  Using  the  same  approach  and  same  unit  vectors 
as  in  (17)  and  (18)  we  can  once  again  set  up  the  following  diagram: 


Figure  21 


A 

From  the  diagram  we  note  that  when  I  ■  S  the  angle  (/T,i)-  4.  the 
right  ascension.  Thus  the  sun  line  and  its  right  ascension  now  may  be  used 
as  a  reference  in  the  representation  of  a, vector  initially  expressed  in  the 
geocentric  system  of  base  vectors  I,  J,  ft  in  the  fixed  system  i,  j,  k. 

A 

Let  M  be  a  unit  vector  initially  expressed  in  the  base  I,  J,  It. 

A  ^ 

M  ■  I  cos<»H  cosl^  ♦  J  cos*^  *inXH  ♦  F  sini|»H 
* 

where  i|»H  is  the  angle  between  M  and  the  equatorial  plane  of  the  earth, 
and  is  the  azimuth  of  M  with  respect  to  the  Greenwich  Meridian  Plane. 
Since  \  and  j  are  in  the  Equatorial  plane  of  the  earth  we  may  write 

I  ■  i  cos  (wt  ♦  x)  ♦  j  sin  (wt  ♦  x) 

J  »-i  sin  (wt  ♦  x)  ♦  j  cos  (wt  ♦  x) 

If  the  time  tm  represents  the  time  at  which  the  Greenwich  Meridian 
Plane  transits  the  sun  line,  then 


i 


•>x 


i  cos4#  ♦  J  sin4# 


i  cos  (wtn  ♦  x)+j  sin  (wtm  ♦  x) 


(23) 


(24) 


b 


! 


♦f  -  wtm 


(25) 


r»‘*  i  •» 
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and  i  =  i  cos [ai(t-tm)  ♦  <t#]  +  j  sin[o>(t-tm)  ♦  *sJ 

(26) 

j  «-i  sin[u(t'tm)  ♦  *  ]  +  j  cos[w(t-tm)  ♦  4>SJ 

which  is  the  same  as  was  shown  in  (17)  and  (18)  .  Substituting  (26) 
in  (23)  ,  'we  get 

* 

M  ■  i  cos^H  cos[w(t-tm)  ♦  *s  +  4  ^  cos<h  »in(ui(t-tm)  +  $s  +  XH1  (27) 

+  k  sini^j 

However  equating  (27)  and  (14)  we  have  ' 
sint^  *  sinG^ 

f 

k 

cosi)^  cos[w(t-tm)  +  ♦  \H]  *  coseH  cos«H 

I 

cosi^  sin[w(t-tm)  ♦  *s  +  *H]  ■  cosOH  sin*(J 

A 

a>^H  “  ®H  ^an8^e  between  equatorial  plane  and  M) 

(28) 

■>*H  ■  u(t-tm)  +  *s  +  ^ 

Now  using  the  program  listed  in  the  appendix,  the  magnetic  field  can 
be  expressed  as 

M  "  ceE  F  +  f  -  \  f  (29) 

where  r^,  Gg,  *E  are  the  Geocentric,  coordinates  of  the  satellite. 

In  the  rotating  system  the  vector  i^is  exp res sable  as 

Rg«ll  cosOg  cos#E  +  j  cosOg  sin#E  +  k  sinGE}rE 
e  is  positive  east  and  is  tangent  to  the  circle  of  latitude, 

L 


•FOUGERE,  P.  private  communication,  L.  G.  Hanscom  Field,  Bedford,  Hass. 


(*■<  I  <* 
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e  ^  is  tangent  to  the  arc  of  the  great  circle  going  through  the  polar  axes 
E 

and  in  the  direction  of  increasing  latitude,  e  is  from  the  center  of  the 

•  rE 

earth  outward. 

3Rp  .  .1  3Ru  1  3Rp 


3  rc 


iSli 

30. 


TgCOSOg  s*E 


In  (29)  the  terms  X,  Y,  Z  and  F  are  given  by  the  referenced  program  and 

F  equals  the  total  field,  X  the  component  of  the  field  in  the  e.  direction, 

°E 

Y  the  component  of  the  field  in  positive  e.  ,  and  Z  the  component  is 

Vj, 

radially  downward. 

Since  the  angles  0£  and  $£  are  the  latitude  and  longitude  respectively 
of  the  satellite  in  the  goecentric  system  (obtained  from  the  ephemeris),  we 
now  have 

er  =  i  cosSg  cos*£  +  j  cosOg  sin$E  *  k  sinOg 
eQ  —I  sin0E  cos*E  -  J  sinOg  sin$E  ♦  ic  cos©E 
e.  =-T  sin*.  +  3  cos#_ 

?E  E  t 

Substituting  (30)  into  (29)  we  find 

A  M 

F  M  *-i[(X  sinOg  ♦  Z  cosOg)  cos*g  +  Y  sin$g] 

-3[(X  sinOg  +  Z  cosGg)  sin*E  -  Y  cose£] 

♦k[X  cosOg  -  Z  sinOg] 

2  2  2  2 
where  F  ■  X*  ♦  Y*  ♦  Z* 

From  (23)  and  (28)  we  have 


1(30) 


(31) 


M  -  i  cos0H  cosXH  ♦  j  cosO^  sinXH  ♦  k  sinOH 
Equating  (31)  and  (32)  we  find 

6^  ■  arcs in  ^  cos9£  ■  Z  sinOgj 


XH  ■  arc tan 


Y  cos<E  -  (X  sinOg  ♦  Z  cosOg)  sintg'j 
-Y  sin*g  -(X  sinOg  ♦  t  cosOg)  costgi 


(32) 


(33) 


/ 


«'•«  I  >* 
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Substituting  (33)  into  (27)  with  ^  *  0^  i.e. 

M  ■  i  cos0H  cos[w(t-tm)  ♦  $s  +  +  •*  COS0H  sin[<‘»(t-tm)  +  ®s  +  *||1 

♦  k  sinG^ 

will  give  the  unit  vector  in  the  direction  of  the  earth  magnetic  field  in 
the  i,  j,  k  fixed  system.  Equating  (34)  and  (14)  we  now  can  find  the 
angles  ©H  and  The  expression  for  angle  0^  is  given  in  (33)  and  « 


=  w(t-tm)  ♦  *  ♦  Xu 

n  S  n 

where  X^  is  given  also  in  (33) .  The  magnetic  field  of  the  earth  is  then 
given  by 


M  =  F  M 


where  F  is  given  by  Fougere's  program  as  a  function  of  (rE,  Og,  4>f) 
F2  ■  X2(rE,  0E,  $E)  +  Y2(rE>  0£,  «£)  ♦  Z2(rE,  0£,  4>£) 


(35) 


> 


/ 


1 
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CHAPTER  V 

DETERMINATION  OF  0  AND  *  FCfl  A  FIXED  SYSTEM 

In  this  chapter  we  will  discuss  the  determination  of  0,  the  angle 
between  the  nose  axis  of  the  satellite  and  the  equatorial  plane  of  the 
earth,  and  $,  the  aximuth  of  the  nose  axis  of  the  satellite  with  respect 
to  the  vernal  equinox.  « 

A.  Theoietical  Description 


The  unit  vector  er  in  (15)  can  also  be  expressed  in  the  form 

A  A 

^  a  uc  m  X  S 

er  "  *  +  BS  +  TmTsT 

2  2  2 

and  dotting  e  with  itself  a  ♦  B  +  T  ♦  2aBM-S  ■  f 
a  »r 

where  M  and  S  are  as  in  (14)  and  (13)  respectively.  If  we  dot  e 

A  A  A  A  r 

with  M,  S,  and  M  x  S  we  find 

A  A  A 

M-e_  *  cosB.,  ■  a  *  BM*S 
r  rl 

A  A  *A 

S*e_  ■  cosB.  «  aM-S  ♦  B 
r  s 

N  X  S»e  - 

|m*s| 

Equation  (39)  can  be  rewritten  as 


(36) 


(37) 

(38) 

(39) 


i 


M  x  S*e 


|mxs| 


The  angle  Be  Is  defined  *n  Chapter  II,  and  is  defined  as 


HZ 

cosBh- 


(40) 


(41) 


i 

i 

i 


i 

i 


fttm 
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where  is  the  component  of  the  earth's  magnetic  field  parallel  to  the 
nose  axis  of  the  satellite  and  H  is  the  total  magnetic  field  determined 
from  the  three  magnetometers  X,  Y,  and  Z. 

i.e. 

H  2  =  X2  +  Y2  ♦  Z2 
0 


(42) 


Solving  (37)  and  (38)  for  a  and  8  by  Cramer's  Rule  we  find 

ft  % 

cosBu  -M-S  cos 8 

H  S 

*  - - r-ri - 

1-(M*S; 


cosB  -  M*S  cos6„ 
s _ M 

l-(M*Sj 


(43) 


a  *  2  •  •  2 

Let  l-(M'Sj  be  replaced  by  |MxS|  ,  then  we  can  express  e^  as 

(cosB..  -  M-S  cosBJM  (cosB  -  M-S  cosBJS  (MSe  )  MxS 

.  _  _ n _ s  _ s _ n  r 

"  A  A  <5  11  '  A  »  ^  “  A  A  1 

T  |||U(*  1^  lu^r  I  !L.r*  I  ^ 


|Mxs|‘ 


|mxs|“ 


|  MxS |  ‘ 


(44) 


Since  e^-k  ■  sin0  where  e^  is  as  in  (15),  then  using  (44)  for  er 
and  taking  the  scalar  product  of  e^  and  k, 

A  A  A  A  A  A 

(cos6u-M*Scos6  JsinOu+CcosB  -M'ScosBu)sin0  +(MSe  )cos0„cosO  sin (*-♦„) 
n  s  n  S  nsrtissn 


sin© 


l-(ft-S) 


(45> 


All  terms  in  (45)  are  now  knownexcept  for  the  triple  scalar  product 

A  A 

for  (MSep.  If  we  take  tfle  scalar  product  of  (44)  with  er  we  get 


(cosBjj  -  M*S  cosBs)  cosBh  ♦  (cos6s  -  M«S  cosB^)  cos6s  +  (MSer) 

1  ■  ■  — - — . "  A1  »  X  .  -.ii..  —i  — — — 

| MxS  |2 


(46) 


wi 
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Replacing  |MxS |  by  l-(MS)  ,  (46)  can  be  rewritten  as 

(MSe  )2  =  l-(M-S)2  -  cos26u  -  cos2fi  ♦  2M-S  cosB,  cos8u  C47) 

I  MS  S  n 

2  2 

Replacing  cos  by  1-sin  8H  and  completing  the  square,  (47)  can  be 
rewritten  as 

(MSe  )2  ■  sinZ8„  -  cos28.  -  (M-S)2  ♦  2M-S  cosB,  cosB„ 

r  n  S  S  tl 

2  2  2  2 
-  cos  Ss  cos  eH  ♦  cos  es  cos  8h 

(MSer)2  ■  s'in2^  -  cos28s  (1-cos8h)  -  (M*S  -  cos8s  cos8H)2 
(MSer)2  =  sin2^  sin2fis  -  (M-S  -  cosBs  cosB^)2 

Now  replacing  (48)  into  (4S)  the  angle  0  can  be  determined,  however 
an  ambiguity  arises  due  to  the  term  | 

(MSep  =*  /sin2BH  sin78s  -  (M-S  -  coss,  cosB^  («9) 

A  A 

According  to  Report  AFCRI, -85-516  page  5,  the  positive  sign  for  (MSe  ) 

A  A  A  A  r 

must  be  taken  where  e  is  on  the  same  side  of  the  M-S  plane  as  MxS,  and 

*  A  A 

the  negative  sign  where  e  is  on  the  opposite  side  of  the  M-S  plane.  Also 

r  a  a 

for  a  flight  where  the  magnetometer  data  is  accurate  (MSer)  will  be  real, 
i.e. 

(MSe  )2  ■  sin28„  sin2S.  -(MS-  cosB.  cosB„)2  >0 
r  ns  s  n 

A  A  2 

In  actual  flight,  however,  it  occurred  that. (MSer)*  <0  which  is 
physically  impossible.  This  was  a  result  of  erroneous  magnetic  field 
data  which  occurred  frequently  during  the  flight  of  0V1-5. 

A  A  A  A 

When  e  makes  an  angle  <90*  with  the  vector  MxS,  then  (MSe)  >0 
and  (MSey)  <0  when  the  angle  made  >90*.  In  the  case  of  a  rocket  flight, 
this  presents  no  problem  but  for  a  satellite  we  are  unable  to  predetermine 

A  A 

the  position  of  ef  relative  to  the  M-S  plane.  In  order  to  get  around  this 
problem,  the  output  for  angles  6  and  6  from  the  program  ASPECT  was  analyzed 
and  those  values  which  gave  the  smoothest  curve  fit  were  selected. 


(48) 


l 


I  m 
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m  ,\ 

In  this  manner  we  determined  whether  to  take ♦(MSep  or- (MSep  during  a 
specific  position  of  the  flight  or  a  combination  of  both.  That  is,  when 
a  crossover  in  the  plot  of  0  occurs  the  sign  of (MSe  }  should  be  examined 
to  insure  a  smooth  fit.  If  at  any  time  in  the  flight  Sj  *  0,  the  ambiguity 
does  not  exist  since  at  this  time 

0  -  0S 

and  #  *  ♦  4 

as  can  be  seen  in  revolution  957. 

Using  expressions  (14)  and  (13)  and  forming  the  scalar  product  of 

each  with  e  we  have 
r 

cos©  cos©H  cos(#-*H)  +  sin©H  sinO  =  cosB(1  (50) 

cos0  cosGs  cos(*-#s)  +  sin©s  sinG  ■  *cosBs  (51) 

Now  according  to  Report  AFCRL-63-871  page  5  and  6,  upon  multiplying  (50)  I 

by  sin3s  and  (51)  by  sin©H  we  can  eliminate  sinO  from  (50)  and  (51).  In 
a  similar  manner  cos0  can  he  eliminated  and  the  above  equations  can  be 
rewritten  as  *  • 


bj  cos©  cos#  +  b2  cos0  sin#  »  a 

(bj  sin0  -c)  cos*  +(b2  sino-  c2)  sin#  ■  0 


where 


a  ■  cosBh  sin0s  -  cosBs  sin0H 

b.  ■  cos0„  sin©.  cos*„  -  cosG,  sin0„  cos#, 

b2  ■  cos©H  sin©s  sin#^  -  cosOj  sinG^  sin#s 

c.  ■  crs0„  cosB,  cos#u  -  cosG,  cosB„  cos#, 

1  It  S  n  S  II  5 

c2  ■  cos©H  cosBs  sin«jj  -  cosOs  cosB^  sin#^ 

the  solution  of  (52)  is 


sin* 


a(bj  sin0-Cj) 

cosOfbjCj-bjCj) 


(52) 


/ 


« 
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or 


cos  * 


tan* 


-a 

cos 


/b2  sinO  -  c2^ 
©  ''blc2'b2cl  ' 


basins  -  Cj 
-Casino  -  c2) 


Another  and  less  tedious  method  of  determining  the  angle  *  is  to  take 
the  scalar  product  of  e^  in  (IS)  and  (44),  both  with  i  and  j.  In  this  case 


tan* 


-  j‘*r 
i»e_ 


where  the  expression  for  e^  in  (54)  is  that  found  in  (44). 
B.  Program  to  Determine  S  and  *  with  Explanations 


(S3) 


(54) 


The  following  Fortran  program  named  ASPECT  was  written  for  the  IBM  7094 
coiqputar  and  its  purpose  is  to  calculate  the  angle  0  and  *  as  described  in 
(45)  and  (53)  respectively.  Frequently  during  the  flight  of  0V1-5,  we  found 
poor  magnetic  field  data.  As  a  result,  when  running  this  program  considera¬ 
tion  of  the  angle  ■  (|M,  e^)  sometimes  produced  erroneous  values.  That  is, 
terms  in  which  cos^{,  written  in  this  program  as  CBETAI1,  were  involved 
produced  at  times  negative  values  for 

'  1 

1.  -  STHETP  * ' STHETP 
1.  -  STHETN  *  STHETN 
used  for  the  determination  of 

CTHETP  -  SQRT(1. -STHETP  *  STHETP)  1 

and  CTHETN  ■  SQRT(i. -STHETN  *  STHETN)  . 

A  A 

In  ASPECT,  CTHETP  ■  COSO  when  the  plus  sign  for  (MSey)  was  used  and  CTHETN  « 

COSO  when  the  minus  sign  for  (MSe  )  was  used.  In  the  final  analysis  these  few 

4  r 

poor  data  points  were  overlooked  to  insure  a  good  curve  fit. 

The  output  of  this  program  is  as  follows: 


**  2  2  A  *  2 
MSe  t  ■  sin  0U  sin  6.  *  (M*S  -  cos0„  cos0„) 

T  MS  5  n 

MSe,  ■  /  MSe  t 
r  r 


>1 


t 


I 


I 


I 


I  I*1*** 


1 


I 
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i  t 

THETAP  and  PH  IP  are  the  angles  0  and  ♦ 
when  ♦  /stSe  t  was  used 

r 

THET^N  and  PHIN  are  the  angles  0  and  *  , 

When  -  /MSe  t  was  used 

r  ;  1  , 

H  =  the  total  magnetic  field  determined  from  the  data 

I  1 

SUNAX  =  (iS,er)  , 

i  «» 

SUNTHE®  (lS,eQ) 

SUNPHI*  (iS.ep 

GMTT  =  Greenwich  Mean  time  of  each  data  ppint 

THETIW  and  PH1NEW  are  the  latitude  and  longitude 
respectively  of  the  satellite  with  respect  to 
Greenwich  (6r  the  time  GMTT.  i  1 

i  i  i 

SIG  is  the  signature  voltage  of  that  sun  sensor 
which  is  giving  sun  data  1 


■  ■  i 

i 

I  ,  } 

I 

'l  '  I 


4 


h 
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tID 

HCP 


MARCOTTL  ASPECT  E3 

IlME=02.PAQLS=i5 -  - - 


ImLL  CONTINUE 

SloJUo  .0 LOGIC  _ I _ 

nerrc  aspcct  list. rep. dfcx.sdd 

C  A5PEC1  - - - 

c  TrlFT  As  ANGLE  RFTWFF.N  AXIS  AN*  EQUATORIAL  PLANE 

C  PHIWIKUTH  OF  AXIS  OF  SA1U.UIJL.  UiBI.  VERNAL-  IQUJNQX 
C  THFTA-MsANGLF  BETWEEN  EQUATORIAL  PLANE  AND  MAGNETIC  FIELD 

X-  PHl-nsAAIKUlH  OF  KAGf.K  1 1 C-  E-IELO  uSJ~.  VtkNAt  LUilUQX. . 

0  IMt.fo ION  TMETAI258I  .PMII258)  .XI 258)  .Y I  258  )  *Z (  258 )  »F  1 258 1  .GM 1 1  ?:>»  I 

UKEG£s2.*3«14l5927/bt>16A,OV4 - - 

9b  READ  ( 5 .90 )  THETAS. PHIS. TM.NuMR.CHECKA.SPAN 

9w  FORMAT  t?F10.S.H0.2.6X.l4.2F6.1  ) - - - 


JO  • 

-  _  .  J*J+I  -ii  .  . . . ,  - 

T«E IAS= THETAS*. 0 174533 

-  PnlS=PHlS*»01JA633  - - 

READ  I  5*2)  I  THE!  AIII  .PHI  ID. XI I  )  »Y  I  I  I  *2  I  II  *F  (  I)  .GMT  I  I  )  .I.l.N'JMQ) 
- 2  FORMAT  1 2FJ-0.3-.1QX  .  5F1-0.3J - 

4  READ  15.11  GMTT.XX.YY.22.SA1 .SA2.SIG.MTEST 

1-FORMAT  -1 7F  lO.l^bXj  15 1 . . . 

C  DETERMINATION  OF  SUNAX.SUNPHI .JUNTHE 

A*ABSISIG-r4.69l  - -  . 

B*ARSI SIG-3.90 I 

C«AB5 1  SI G-3. 1 9 )  - - - -  -- 

D=AbblSIG-2.4~ )  ■ 

£=Abb  (SIGMAS) .  - - 

FF=A6SI  MG-1  .57) 

IF!  .31-A4 — 5-*10.10  - 

5  IF  (.35-8)  6.11*11 

6  IF  1.35-0  7.12*12  - 

7  IFI.40-DI  8.13*13 

-  8  IF  l.3S=£FA_9.44.14 - 

9  IFI.35-EI  97.15*15 

lu  ALPHA1 =-45. 1 13-1 .  27  3*SAl  +  l<U.39.9.fi5A-l.f.SA.l  -2-»lS7*-SAl--*3 —  . 

ALPHA? =-44. 852+1 . 363*5A2+ 13.063#$A2*SA2-1.984*SA2**3 

_  SUNPHI s90.-ALPHAl  - 

SUN THE *90. ♦ ALPHA 2 

- ALP11A1  -ALfiHAl  *-01 74533 _ 

ALPHA?* ALPHA 2* .0 1 74533 

XCOb=SOKJ  (COS  I  ALPHA  I )  »*2-SIN(  Al.  PHA2  1**2  I _ _  .  . — 

XSIN=S0RTI1.-XC0S*XC05) 

SUNAX=ATAN(X5IN/XC0bl . .  . . . 

SUNAX»SUNAX*57. 29578 

_ GO.  JO  46 _ -  — - _ _ _ 

11  BETA1 =-45.956  +  4. 636*SA  1  +  1 1.970*SA1*SA1-1.901*SA1**3 

.  -  .  BETA2  =  -45. 291  +  3.  724»SA2tl2.0.UtSA25SA?-I  .a73ttSAZ*Jfc3-  : _  ....  _  — 

SUNPHI=90.-RETAl 

...  SUNTHE=50.-DETA2  - - - 

BETAloETAl*. 0174533 

- bELAZ=f>£IA2*«01  74533  - 

XCOS  =  -SOKT(COS(PETA1  )«*2-SINIRETA2)**2) 

- .-XSIN=bQRTin-XCGS*XCOS) - - 

SUN AX  =ATANIXSIN/XCOS) 

SUNAX=SUNAX*57. 29578  +  180, - ! -  .  . 

GO  TO  16 

— L2-.GAMMA  =  -45»8-7-7*2, 53TV*SA  ]  4 1  fc.6Xfet5Al.lS  Al  r  1 . 994*  5A  1 1»3 - 

SUNAX=90. -GAMMA 

SUNTHfcs90.-GAV.MA2 


6AVMAxC'AMyA».0174S13 

GAMMA2=GAKMA2*.01 74533  - 

XCUS=  SOK  T  (  COS  ( GAVf-  A  )  **2-SI  N (  GAMKA2  I  **2  ) 

-_XiINs^aaT4l.^XO!j»XCOi>  - 

SUNPHI=AIAN(xSIN/XCOS) 

5UNPHI sSUNPHI X57.29578  - - - 

GO  TO  16 

13  F>FL  T  A=-45.  795*3.  296»SA  1  ♦  12. V26 *SA1  +&A1» 1  ^970‘SAl  **3— 
SUNAX=40. -DELTA 

0ElTA2=-45.5g4*3.66ti*  jA2*12.1CT.»SA2lSA2c4^m-Sfc»A2**2  . 
SIJM  rC=9P. -DELTA? 

DEL  TA=0LLTAJU01  74533  -  -I - 

DEL TA?rOELTA2*. 01 74533 

XC0S  =  -S'JR1  (C0S(JELTA)«*2-sH±U)El.JA24**24  — - 

XS1.\  =  SGRT(1.-XC05*XC05| 

SUNPH 1  sATAN (XS1N/XCQS)  *57.29578+180.  . 

GO  TO  16 

14  FTHFTAr-45.6n]+2.5?3«SAl  +  l?,483**A4*SAl«4.Q}<M&AV**-3- 
SUNAX.-9fi.-F  THETA 

FTHET2  =  -46.021*3.?28#SA2*12.M7fcSA2*»2-  l  .A95*SA2*a3— 
SUNPH 1  =  90. -F THE  12 

FThLJAsFIhETA*. 0174533  _ 

FTMCT?=FThET2*. 01 74533 

XCOS=-SUKT  (COS  (F  THE  IA)**2-SIN(FTMFT2  1**2  l  ........ _ _ 

xsiN=soRr(i.-xcos«xcos> 

SUN  1  Ht  =A1 AN( XS IN/XCOS)*  *7.295  7-8+  1FO.  - - 

GO  TO  16 

15  EPSlLNs-4e.057+l4.876»SAUS.»l5Al.SAlJt.SAl.^.&&4»SAU.«.3  ... 
SUN Ax  =  90 .-EPS  I IN 

EPSIL2=46.387-7.4?4*SA.~!-9.  191*SA2*ft2*1.4kOO*SA2.#^ - 

SUNPH!r90.-EPSIL2 

EPSILN=EPSILN*. 0174533  ...  _ 

EPSIL?sEPSIL?*.01745,3 

XC0S=S0RTIC0S(EPSILN)**2--SlM-EB£.U24.t*24 - 

XS I N=SQRT ( 1 .-XCOS*XCOS1 

SUrN THE.*-ATAN(XSIf«/XCUS)*57.29578  _  _ _ 

DETERMINATION  OF  BETA-H 

16  XMG=249.45B*XX-604. 909  _ 

YMG  =  251.2 56 *YY -606. 030 

ZMG=247.934*ZZ-592. 562  _ : _ 

H=SGRT ( XMG*XMG+,YMG*YMG+ZMG*ZMG) 

GO  TO  30  _ _ 

26  J  =  J+1 

30  REFER=GMTT-GMT(J»  _ _ _ • _ _ 

IF  (REFER)  23.24.24 

23  REFER  =  -REFER  ....  .* _ 

24  IF  (REFER-CHECKA)  25.25.26 

25  T I ME=GMT  T-GMT  (  J)  _ _ , _ 

TIMET=TIME/SPAN 

IF  (TIME)  27.28.29  •  . . . . 

27  FNEW  =  F(  J-J-TIMET*  (F(  J-l  )-F(J)  » 

XNEW  =  X( J) -TIMET* (X(  J-l  )-X<  J.U — 1 - 

YNFWrvj J)-TIMFT*(Y(  J-l )-Y( J)  ) 

ZNEW  =  Z(  J)-TIMET*(Z(  J-l  )-Z  LJj-L - 

THE  TNw'r  THETA  ( J)-T  IMET*  (  THE  T  A  (  J-l  ) -THETA  (J)  ) 

PH  I  Nt/;  =  Ph  1 1 J  )-T  I  ME  T  *  ( PH  I 1  J- U-^P-M-LU-U _ 

GO  10  33 

28  FNEW  =  ftJ)  _ _ 

XNEw*X ( J ) 

YNEw  =  Y  ( J » ..  _ _ 

ZNFW=Z( J) 


a 


TrP*'Jnlsrr<E  T  A  (  J  | 

i  JV  - 

v*..  rr.  33 

Sfi.-  *:f  t  Jt-r  J.)  — r  l  j*n  I, 

»*  -  «=x  j  j  »-r  j  vc  r- : » :  j  i -x  i  j+i  i  > 
v\fMsY|^).-ri»ri':  vtj.t--Yt.jAtw- 
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£\-:  .v=i  c  J  J  -  ” !  v£  r  ;  c  I  c  J  i  -  £  c  j  *  1 )  t 

Ir*E  l^Mz-1  ;n«i  Z.A  t-JI  -  I  UtlXtl 

^-!i:ic.Ysr-t-<i  i  j  i  -  r  :•••;  r*i a-i  c  j)-phi  i  j-*-i » » 

A3  Cai  LAitiXAGi/H  _  .  - - - 

T ►»£  r!\(/riT rl£  TN  a  1  «  ‘  3 5 

Ph  IN£.i'=8n.iMiA 1 1 ’Aiii - .  .  - - 

Ct  *£f'vlNAT  Iw>i  r--  T A-H  AND  PhI-h 

.1  T nr  7h=W Xf  *. > r .V >  l- 1 Kt  Ux.pUil{4i&tk!  It*r  t.VyJ  I  X FAl*  4 - ... 

C"n£  fHa  jiJRT  (  1  »-C  TeF  T h-*;>TH£ TH ) 

IhF  TAK=A.IAAi.(  iklHE  Ir./CIttE-lHJ - - 

r  ;  m  s  v  N  £  s  *  i  i  r » •  f  *  \ \  r  •*  *  C  C  5.  (  7  r-  F  7  '  <  y  I 

AM  J  Aras  A  T  «u\.I 

IS  INC  Pm INFwt -TER* **:C><?mINE*>  »  ) 

vn ;  HsCjiESAJurv-.it- t.v. )  lA^Ajiraa - - - 

OF  TpI^-iINATIi.N  Of  ANO  MSER 

.  Adf.=  COST. THETAS. I  _ 

INC  THFrAS) 

ArP=CCS(.PHlo>  - -  - 

OFHablNCPHir.l 

FMOxU.VrC  TrfF  7H*Ar.C  *C.^S-tPh4Ji»PHUAW^yu«F4Ji»£FF - 

G'iMA>  =  £  JNAX*..*  17453? 

ElNAJj^U^lStUiUUUJ - - 

■I'.:  iUNsfOSCGuNAX  ) 

Cil  jBf  T=CGJAJI»CF.£  t  Ah  - - - - 

eMir*T  =  SINSUN*GIN'XN*(  l.-Cf-ETAH*CPETAH!  -(E*DOTS-CO'RET  )■ 

IF  {.EM.SEi.tt-  6-1  *-6 -1  *6-2 - - - 

EviERsO.O 

-00.  tO-4J - 


oi 
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uMsCPsSwWT C  £Mj£R T ! 

DETEST*  I.  NAXiON.  CF  THETA  AND  P-fcU- 


AA=(OEF-EMDCT6*orHETH)»COSUN-*-CSTHETH-EMDOTS*DEF)*C5ETAH 
B*i^AiaC*ClME.tH.*SlN.  l.PH.j. k — - - -  . 

CC* 1 •-EMDOTS*EMCCTo 


SThETN=(AA-RR)/CC 
CTHE-IPsEOR-T.  C.  U--S.7.HETP  *S.tHEXP4.- 


THE T AP=AT AN (SThETP/CTmETP 1*57. 29578 
-  ClHETMsSas.I  CL«.-S.lHEtN.»E.lM£INt - 


THETAN  =  ATAN(SThr.TN/CTHETN)*5  7.29578 

- E£tUt u  CQS.C. £  tiUtl  .  -  *  - - 

SPHIHsSINJPHIH) 

fcll*ClHETH*Q£f  *CPnIh-iMlC*.EIn£JLtltAC£___ 
R2=CThfTH*DEF*SPhIh-apc*.SThFTH»DFH 
Cl=CItiFTH*CQSUN»CPHlH-A&C*E£EXAJ±!tAjC£_ 
C2sCTHFTH*C0SUN*SPh:h-APC*CRETAH*DFH 

..  .  AMilM.P.t»STHFTP-f  1  . . . 

AI5ENOM«-(  R2*STHETp-C2  ) 
BNJM=bl>wTrl£TN-Cl  .  _ _ 


72 

73 


PCFNOMr-( R2*5ThE  7N-C2  ) 

PH  I P = AIAN.  1  AjNUM/ACENOK  l  »67.A967fl 
H  -t  I  N  -  A  T  AN  t  n  NO/  aGE  No-*' )  *5  7  •  2  95  7  8 

IF  ( PuIELl  JZQ»  72.*  72  _ 

IF1AMLMI  73*75*75 

p‘-!ip=PHi?.tj.ac.  _ 

GO  ro  76 


Ifj 

76 

82 


7v,  I F  (A\lV )  71.73.73 
71  PhII>*Ph!P*360* 

GO  TO  76 

PrilP-PHiP _ 

IrlPHIM  80.82.82 
I F  ( bhuy. )  63  .83*86 
bi  PmIN.-PhIN*180. 

GO  TO  44-  ------ 

tiu  IF  (bNUMI  81.63.83 
81  PHIN~PHlfc*866«— -  - 
GO  TC  44 
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bi  PHlbsPMltt  - 


44  THETNVt  =  1rtFTNW*57.?Ob7fl 

...  PHJ«EK=PHIN£W*3Z.2a678 - 

taRITI  <6.49!  EMSEM .EMSERT *H. THE TAP. THE T AN, PH I P.PHIN .SONAX.bUNTn 

.  .2suaam-.  Inman. pm.1 sitAtSi u>. uaj.* - 

49  FORMAT! 1X.2F7. 4. lOF 6.1 .F6.3.F  10.31 


r 


lPMINEw.GMTT 


-_43  FORMAT! 10F6*4 ,F10. 31  — 
97  IF  ! MTCbT-999 I  4,96.96 

_ S6.  CALX.  EXIT _ _ _ 

STOP 

END 


SDATA 


■i 
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CHAPTER  VI 

DETERMINATION  OF  (fe^,  U") 

In  this  chapter  we  will  determine  the  angle  between  the  unit  vector  U" 
as  expressed  in  (19)  and  e^  in  (22).  The  unit  vector  e^  is  the  direction 
of  sensor  C  on  the  satellite  and  e^  «-e  x  ep  for  the  e  ,  eQ,  e^  systtn  • 
with  regards  to  the  satellite  as  discussed  in  Chapter  II. 

A.  Determination  of  e,  in  another  fixed  system 

The  problem  as  was  mentioned  in  Chapter  III  is  to  solve  for  the 
angle  P.  In  order  to  do  this,  we  must  set  up  another  system  for  e#.  In 
the  following  figure,  e+  is  in  the  plane  cf  ^  and  N2- 


1er 


Figure  22 


and  can  be  expressed  in  the  form 


e.»  aS  ♦  ®^er  X.S^  ♦  x  S) 

♦  ,  TT-rfir  TS^eV*S)T 


where  S,  and  the  vectors  (er  x  S)  ^  S  x  (er  x  SJ 


|er  x  S| 


IS  x  (er  x  S)| 


define  an  orthogonal  system  of  unit  vectors.  Since 

S  x  (ef  x  S)  ■  ey  (S-S)  -  S(er*S)  *  er  -  Scos6s 
it  is  easy  to  show  that 


(SS) 


rfl 

i. 

li 

•’h 


f" 
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|cr  x  S|  *  |S  x  (cy  x  S)|  *  sinf»5  (S6) 


where  t$s  «  (fS  e^) 


To  determine  a,  |5,  and  y  in  (55)  we  must  take  the 

scalar  product  of 

e.  with  S,  e_  x  S  ,  and 
$  r 

S  x  (er  x  S) 

respectively.  This  gives  us: 

sinBs 

slnBs 

a 

a  »  e..S  »  cosy. 

^  5 

(57) 

«  V  <er *  *’  . 

Ce4«rS) 

(58) 

sin6# 

sinBs 

y  *  e..  ^er  x  «) 

.  *V  i!l_ 

-  S  co$6  1 
s 

C0SYS  COS0S 

,  (59) 

*  sinB 

sir  6. 

sinB, 

f 

s 

s 

s 

1 

4 

Using  (57)  ,  (58)  ,  (59)  it  is 

now  possible  to  set 

• 

m 

• 

•  A 

e.  ■  cosy.  S  ♦  ^e4er^ 

(®r  *  S)  - 

cosgs  cos^s 

S  x  (er  x  S) 

(60) 

*  5  sTOs 

sinBs 

sines 

sinBs 

If  we  consider  now  the 

scalar  product  of  e^  ir. 

(CO)  with  itself. 

1  -  COS2Ys  ♦  <heT  S)2  ♦  COS2Y  COS2Bs 

sin20  .  sin2B 

s  s 


>: 


t 


so 


2  2  2  2  2  *  2 
■  >  sin  8S  ■  sin  8S  cos  ys  ♦  cos  8S  cos  ys  ♦  (%er  S) 

*>  (e^e^S)  *  *  i4in^8s  -  cos^y^ 


(61) 


Inserting  (61)  into  (60)  ,  e  is  now  expressible  in  the  orthogonal 

A  A  A  »  A  A 

systea  S,  er  x  S  ,  S  x  (er  x  S)  .  The  plus  sign  for  (e#ey  S)  will 
sin8  sin 8 

9  9 

A  A  4 

used  when  e.  is  on  the  same  side  of  the  e  -  S  plane  as  e_  x  S.  Otherwise 
the  sinus  sign  will  be  used  in  (61)  . 

B.  Determination  of  the  angle  p 

Using  expression  (22),  (60)  with  the  value  of  (e6e'rS)  found  in  (61 ), 
the  purpose  of  this  section  will  be  to  deteraine  the  angle  p  that  e.  rakes 

A  A  A  ▼ 

with  Nj  while  rotating  in  the  Nj  -  S2  plane.  In  this  method,  we  shall 
equate  the  h  coefficients  in  each  of  the  above  aentioned  expressions  for 
eA.  The  k  coefficient  of  (22)  it  cosG  cosp.  If  we  let 

7 


A  «  cosG  sin«  sin©s  -•  cosOr  sin*s  sinO 

B  »  cosO  cost  sinG  -  cosG  cose  sin©  (62) 

9  S  9 

C  -  cosG  cos©  sin(e  -e) 

5  5 

A 

It  can  be  shown  that  the  k  coefficient  of  (60)  where  S  and  er  are  in  the 
i ,  j ,  k  systems  is 

A 

cosy.  sin©_  *  c(*$erS)  -  C0SYs  cos8s  f'cosG-  cose.  B-cos0_  sine.A)  r6.Vl 
s  s  2  2  •  l  s  5  s  -  I 

sip  6S  sin* 8, 

Equating  cosq  cosp  to  (63)  and  simplifying 
2  * 

cosp  •  sin  8  cosy  sinO  ♦  C(e.e  S) -cosy  cos8  cosG  [  cosO  sinG-  cosGsinO  cos(e  -e) ]  (64) 

_ is  5 _ 5 _ ™  a  5  5  5  9  _  5 _ 3 _ 

2  1 
sin  8S  cosG 


i 
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We  nay  further  simplify  (64)  by  using  the  expressions  for  ey  and  S 
in  (15)  and  (13)  respectively. 

■>e  *S  ■  cose  cosO  cos(t-« )  ♦  sine  sine  »  cos  5  (65) 

r  s  5  5  s 

■>sin^?s 


*  sin*6s 


cosp«  cosys  [sinOs  -  cosBs  sine]  ♦  C(e^er§)  (66) 

sin^Bs  cose 

Due  to  the  tern  (e^efS)  an  ambiguity  results  for  the  angle  p.  As  ' 

can  be  seen  by  figure  IX  in  this  chapter  the  angle  y  is  a  minimum  when 
ef,  S,  and  e^  all  lie  in  the  same  plane.  We  know  from  a  previous 
discussion  that  (e  e  S)  is  positive  when  e  is  on  the  same  side  of 

*  A  ▼  A  ▼  A 

the  er-S  as  exS.  Therefore  we  may  also  say  that  (e^erS)  is  negative 
as  ys  goes  from  its  max  value  to  its  min  value. 

Ne  can  obtain  another  expression  for  the  angle  p  by  considering 

A 

the  scalar  product  of  (22)  with  the  sun  vector  S  expressed  ir.  (13). 

If 


.2 


cos 9  cosp  «  sin  6S  sin©s  ♦C(e^erS)  -‘cosBs(sinO  -  sinO^osBjl 


cosy, 


cosy. 


sin9.  ♦  sin0,  cos*B.  ♦  C(e4e  S)  -  cosS  sin0 
5  s  s  v  r  s 


COSY. 


Aj  *  cosO  sinOs  -  cos0s  sinO  cos(*-$s) 
A2  ■  cos0s  ^in(6-0s) 


b 


Then 
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but 


e.»S  ■  cosy.  *  A.  cosp  -A,  sino 

2  2  2  2  2 

r>Aj  cos  o  -  ZA^sinp  coso  ♦  A^  sin  p  «  cos  ys 

■>Aj2  ♦  A £  tan2p  -  ZAjAjtanp  -  cos2ys  (l*tan2p)  «  0 

*>(k^  -  cos2y#)  tan2p  -  ZAjAjtanp  ♦  A^2  -  cos2ys  ■  0 

— . .  .  . . .  — 

■>tanp  ■  AjAj1  /  Aj  cos  ys  ♦  cos  y,(Aj  -  cos  yj 

2  2 

kJ  -  cosS. 


/a  2  ♦  A  2  - 


t»np  *  *  C0SYS  f  Ajfc  ♦  A2 

T 


cos  Y« 


Aj  -  cos  y, 

2  2  2  2 

Aj  ♦  A-j  ■  cos  0  sin  ©s  -  2 cos©  sinGs  cos©s  sinO  cos(«-$s) 
♦cos2©  sin2©  cos2(e-e  )  ♦  cos20  sin2(*-$  ) 

9  9  9  9 

2  2 

*  (1-sin  0)  sin  0-2  sir.O  cosG  sin0_  cosO  cos(t-t) 

9  9  9  9 

♦  cos20.  (l-cos20) ’cos2(*-*c )  ♦  cos20.  sin2(6-6  ) 

9  9  9  9 

*  -  sin20  sin20  -  2sin0  6os0  sin0  cos0  «»(♦-♦) 

9  9  9  9 

-  cos2©  cos20  cos2(»-*  )  ♦  1 

9  9  • 


■1-  (sin©  sin*s  ♦  cos©  cosOs  cos(6-*s)^ 

so  by  (65),  Aj2  ♦  A 2  ■  l-cos2Bs  *  sin20s 

Inserting  (68)  into  (67)  we  now  have  another  expression  for  angle 
p,  that  is 


- - - - — 

tanp  ■  AjAj  *  cosys  /  sin  09  -  cos  ys 
- 5 - 5 - 

k2*  -  cos  y9 


(67) 


(68) 


(69) 


.  ! 

i< 


7  ' 


( 


S3 


The  ambiguity  can  be  resolved  in  the  sane  manner  as  previously  discussed 
since  the  term  in  the  square  root  is  exactly  that  found  in  (61).* 

C.  Theoretical  description  for  (je»,  U*') 

Rewriting  the  expressions  for  U"  in  (19)  and  e^  in  (22),  we  have 

0"  *  i  cosQg  cos{u(t-tm)  ♦  *#  ♦  *£}  *  j  cosGg  sin(w(t-tm)  *  *s  ♦  *g) 

♦  k  sinO£  •  * 

e^  *-i(sind  cos*  cosp  ♦  sin*  sinp)  -  j(sin@  sin*  cosp  -  cos*  sinp) 

♦  k  cosO  cosp 

Now  taking  the  scalar  product  of  (19)  and  (22), 

*  •  • 
e#.U"  ■-  sinO  cosp  cosCg  cos(*-[w(t-tm)  ♦  *s  ♦ 

-  sinp  cosOE  sin(*-[u(t-tm)  ♦  *#  ♦  *£]) 

*  * 

e#.  U"  -  cos(«e^,  U)  ■  cosGgfcosp  (cos0  taneE  - 

-sinp  sin(*- (w(t-ta)  ♦  *s  ♦  *EJ)J 

In  summary  (70)  determines  the  aag'e  between 
angle  pis  determined  by  (69)  or  (66).  When  yg  is 
♦  sign  for  (e^e^S)  and  when  ys  is  decreasing,  use 

D.  Program  to  determine  (*ea,  U") 

In  the  following  program  written  for  the  IBM  7094  computer  and  named 
ASPECT  FINAL,  the  output  is  as  follows: 


♦E]>  | 

a 

♦  cosO  cosp  sin@E 

sin©  cos(*-[u(t-tm)  ♦  *s  ♦  *E])) 

(70) 

e#  and  U"  where  the 
increasing,  use  the 
the  -  sign  for  («$*rS) 


*  The  expression  for  sinp  and  its  derivation  can  be  found  in  Appendix  F. 
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/ — 2 - T 

ROMP  is  the  angle  p  when  *  /  sin  -  cos  ys  was  used 

/ - 2 - T 

ROWN  is  the  angle  p  when  -  i  sin  8^  -  cos  ys  was  used 

A 

ANGLEP  is  the  angle  between  e#  and  U"  whon  ROW'P  was  used,  and 

A 

ANGLE!!  is  the  angle  between  e^  and  U"  when  ROWN  was  used 
BANGLP  ■  360  -  ANGLEP 

BANGLN  -  360  -  ANGLEN  , 

t 

SUNPHI  ■  ys  *  the  angle  between  the  sun  and  e# 

GMT  ■  Greenwich  mean  time  for  each  data  point 

The  only  problem  that  occurred  while  running  this  program  was  at  varying 
times  the  term 

A 

oin20s  -  cos2ys  |  (71) 

was  found  to  be  negative  and  therefore  an  error  was  encountered  when 

A 

determining  (e#erS).  However  this  occurred  so  infrequently  that  the 
data  points  at  which  (71)  was  negative  were  just  overlooked. 


* 


b 


i 


r- 


*•••«  i 


( 


TIP-'"'  n 2  3?  ”Var£()TTT"  A'r‘:r>TT  FINAL  H 

MCP  TIVC=0?.PAGES*lb 

"'  "JALL  ~  '  CONTINUE  . 

DLOGIf 

Trrr rnrwnn  l i m  .ref  ,i>ecn ,sp73  • 

C  ASPECT  FINAL 

“  C  n*FT  r’R'M  INAT  ION  ‘OF  ‘  Wf.'W  fOI<  ANGLE  BETWEEN  E-Ph'I  AND  U  IN  FIXLD  SYSTEM 

C  IF  SUNPHI  IS  INCREASING  FROM  0  TO  180  DEGREES*  USE  ROWP 

C*  '  IF  SUNPHI  IS  n.CRIASING  FROM  180  TO  0  DEGREES*  USE  ROWN 

_ 9  READ ( 5  *  5  I  THC T  AS.PHI 5*TM _ 

5  FORMAT (2F10.r.F)n.?( 

TMETAS=THETAS*.0) 74533 
PHIS=PHIS*. 0174533 

UMEfa A==2.«3.  1415V 2 7 /86164.Q91 _ 

iT'RTAol  5,1)  THE  T  A  ,PH  I  ,SUNAX,St'NPHl  ,  TFlE  TN.v  ,  PHI  NEW  *GMT  »MTE5T 
1  FORMAT  (6F7.1  ,F)0.?»4X,I4) 

- run  a  =thftkv.  rrwrs - ; - 

PHI=PHI*. 0174533 
S  U  NAX  =  SUNAX*.ni74533 
SUN PH  I sSUNPH I*«0174533 

~A  =  C (TS  (  T HE T  A  )  * S  I N  I  THE  T  AS  ) -COS  (  THF TA S  )  *S I  N  1  THE T A  )  *COS  [  PHI  -PHI S  ) 

RsCGSl THETA S)*S IN (PHI-PHIS) 

- CTCI/STSUNPHT'I - - - - - 

TERM=C*aURT ( S I N ( SUNAX >**2-C*C ) 

SKOV/PsA^B+TTRtf  : 

SROWNsA*B-TERH 

- ClTOv:=h*tI-fVC - - T_ 

HYPOTPrSORT(SROwP*t2+CROVUH»2 )  , 

- HYPO  I  NrGTJRl  (  E!<LmN*#2K«UW*«2  1 - f - 

SRCtvPsSROwP/HYPOTP  * 

- -  ■  SWTCJJ  vrmwv  7177  PTj  t  ^ - 


^  CROi-.'PaCROW/HYPOTP  _ 

c  determination  of  ANGLE  FFTWCEu  e-phi  and  u 
TWA s’ PT7T^Fw*ToT74TT3  " 

PSIeTHCTNW*. 0174533 

■ - Tr*T0 STTFIFTirnr5 TK t  PS  I  I7CUE  ( PM  I -S  ITf  I  THE  I  A  I  *C0S {PHI -IOMEGA* - 

2<GMT-TM)+PHIS+AM0A)I 

. . r «STN"(  PHI -j OMEGA#  (GN'T'-Tf'1 )  *PHI5+AM&A  j1) - 

XCOSP  =  COS(PSI  )*(CROwP*D-SROV.'P*E) 

- jTS T7IFJ5CI? 1 1  i  ;-kG05P#xCDSP  I - : 

ANGLEPr ATAN IXSINP/XCOSP 1*57, 29578 

- Tc^sHKO?TI^"#'Ii:,R^’AN_«6-,sliiO'--N#E) - ; 

XSINNsSORT ( 1,-XC0SN*XC0FN) 

- AWG  LTff=  ATA  RTX3T  N  \  /  X  L  U  S  M  I  #  b  / .  i»  y  b  7  8 - 1 - , 

RANGLPsATAN(-XSINP/XCCSP)*57.29578 

- mTn^rnrTTTS - — 4 — : - 

4  ANGLEP=ANGLEP+1B0. 

-  "  B7ITICO 1  =  n ANGL P 1  W>  • - - 

00  TO  12  U 

- ft  TO  LT«  F  TOGtTTTTOT; - 

12  BANGLNr ATAN I-XSINN/XC0SN)*57,29578 

- I'MXLUSNI  /*  I,b - 

7  ANC,LEN  =  ANGLEN+180. 

- 3  A  N  (3  L  Nr  B  fl'NG'E.  N  *  1  HO . - 

GO  TO  13 

- B"miNGLN»H7rNGLNV3fc'0. - 

13  ROWPs ATAN( SROWP/CROwP 1*57,29578 

- mr  M  *  ATT'A.^'1 5WTOgfflTT?  57729  b  7  8 - . 

IF  (CROWP)  24*3*26 

- Zg"'T7nyp3RfpRfprr?g; - 

I 

_ _ _ • _ „ _ i 


¥ 


b 


l 


i 


_ _ _ _ 56  _ 

'  GO  TO  46 

lb  IK  ( SROWP )  44*44*46 

'  WTtOfcPsR  0  w  KTT60  •  . . . 

46  IF  (CROwN)  *4,3.66 _ 

*4  ROW^  =  K'0.*,N  +  1  8 0 • 

GO  TO  66 

56“ If  (SROwN)  64,64,66 
64  ROWN=ROwN+360» 


i 


57 
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A 

E.  Plots  of  the  Angle  between  and  U"  with  explanations 

When  plotting  (Je^,  U")  the  rule  to  choose  ANGLEP  when  ys  was 
increasing  and  ANGLEN  when  y$  vas  decreasing  could  not  be  strictly 
adhered  to.  The  angle  y  could  start  increasing  or  decreasing 

5  A 

without  switching  to  the  other  side  of  the  e^-S  plane.  In  order  to 
determine  if  the  switchover  actually  occurred,  the  appropriate  angle  p 
had  to  be  examined (ROWP  for  ANGLEP  and  ROWN  for  ANGLEN).  A  switchover 
from  ANGLEP  to  ANGLEN  when  yg  starts  decreasing  should  incur  a  smooth* 
change  from  ROWP  to  ROWN  and  similarly  when  y  starts  increasing. 

*  A 

Another  point  to  bear  in  mind  is  that  unless  (IS,  er)  is  quite 
small  when  a  critical  point  occurs  for  y  ,  then  to  insure  a  switchover 

A  5 

from  one  side  of  the  er~S  plane  to  the  other  side  the  angle  yg  should 
have  a  min  value  fairly  close  to  0*  and  a  max  value  fairly  close  to 
180*.  This  can  best  be  seen  in  figure  #22.  The  vector  e*  lies  in  the 

a  a  •  a  * 

N,  -M?  plane  and  the  min  value  of  y$  occurs  when  er,  S,  and  e^  all  lie- 
in  the  same  plane.  { 

Still  another  point  to  bear  in  mind  is  that  high  detector  readings 
may  occur  on  detectors  looking  in  the  direction  of  sun  sensor  D  even 

A 

though  (4e#,  U")  is  a  large  angle.  This  occurred  in  revolution  480 
approximately  53.5k  seconds  GMT  and  the  reason  was  that  the  detector 
was  looking  almost  directly  into  the  path  of  the  sun  as  shown  by  the 
plot  of  (#S,e^)  for  revolution  480*  If  the  detectors  in  the  direction 
of  the  sun  sensor  D  are  not  looking  into  the  path  of  the  sun  and  if 
there  are  no  reflections  from  the  albedo  as  may  have  occurred  in  revo- 

A 

lution  957,  then  these  detectors  should  have  high  readings  when  (fe^,  U") 
has  low  angular  values  and  low  readings  when  (fef,  U")  has  high  angular 
values. 

As  one  can  readily  sec,  in  the  following  figures  #23-#26,  there 
are  points  which  do  not  follow  the  general  trend  of  the  curves.  These 
stray  points  were  included  in  the  plots  to  give-  a  complete  picture  of 
the  data  analyzed.  The  data  listings  for  these  plots  can  be  found  in 
the  appendix. 
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DEVOLUTION  480 
'J'iT.  ANGLFP  FROM  *851'?  SEC  •.  TC  53414  SECS. 
‘t<f;T7TT  AND  ROwP  1 N  D I C  *•  T  T  SWITCH  l6  THE  ( 
/NOLEN  THROUGH  55896  SECS. 


A'iGLEP 


bANGLP 


A  0.320 


46.462  1 C  6. 1 74  05  3.826 


1-8.82? 
1677077 
2-4.4  16 
2TC7676' 
136.4  1'*' 


97.880 

TT6.T3T 

143.277 

rvrrrrv 

lr.i.in.i 

~7TI".TR2 

1  02,008 


162.20V 
— TTcmrr 
97.966 

- FTiZr~r~ 

U.',T>  /. 
— TT'- 

— rnrvTTo  ■ 


197.791 
■'  2  V  tb  T'i' ' 

262.644 

772.318 

217,246 

TOVTbTB 

207.03P 


i  *. 


[  |  I  I  [i—  L  1|  I  II 


314.082 111.000  466 24.000 
313. 65C  107.900  48655. COO 


5  5. 798  304.202 

61.099  293.901 

102.601  4  r  72  3  .  ■'0  0 

96. 2 00  48:60.000 

62.110  297.890 

63.059 - rrSTSSI — 

92.600  .48  791 .000 

62.147  297,853 
61.101  298.899 
59.608  300.392 

57.245 3o2.715 

87.500  4666*. 000 

84.500  46657.000 
80.20.1  4692S.O0C 
75.900  46565.000 

54.547  305.453 

5TT.T54 - 559.044  1 

74.300  46996.000 
71.20'"'  4  9  0  1 

46.129  313.871 

r*  "ror - rrs-fA'A — 

•  68.500  49065.000 

7V.4"i4  280.596 

75.877  284.123 

737767 - 2‘66.?5?"~ 

72.288  287.712 

37.000  51953.000 
36.300  61964. COO 

’  33.300  'mTT^orT 

3?.8on  52032.000 

1  16.513 

105.655 

157.338 

202.662 

163.360 

196.640 

167.300 

O4* 

5  3  ',46 . 000 

m.?ci6 

105.115 

158.714 

201. 2e6 

161.897 

198.103 

170.500 

53383.000 

1-7.  bi<i 

107.361 

1  58.975 

201.025 

"TTa.  992 

201.008 

175.000 

53414.000 

1-8. 479 

104.477 

157,796 

202.204 

158.169 

201.831 

177.100 

5345*. 000 

■fMfOlfl 

mmrmm 

■mw 

204,409" 

204.^12 

1  7  b  •  4  C>  \j 

33461.000 

1 12.035 

98.911 

155.540 

204.460 

154.686 

205.314 

172.7C"' 

53616,000 

1  i 7,4  /I 

'177593 

T  5  4. 3 88 

205.6  14 

150.310 

209.690 

166.70  ' 

53549.000 

117.663 

93.181 

166.055 

204.945 

149,768 

210.232 

165.600 

53686.000 

3-4.57? 

89.602 

26.034 

331.966 

147.581 

212.419 

162.100 

5  36  17.  COO 

3-2.413 

83.705 

30.438 

729.562 

146.927 

213.077 

158.300 

5  3o64 .000 

5-7.364 

36.09T- 

25.708 

334.292 

141.730 

218.270 

15  5.60.' 

53685.000 

’14.390 

81.545 

27.68? 

332.311 

136.320 

'223.680 

149.50' 

53722.000 

316,770 

TTTTor" 

30,  oo<5 

3  2  9 , 0  9 1 

1 3  4 .  ?  •>  3 

228,602 

146.90O 

33743.000 

2  s  5 .  <1 6  1 

55.56H 

62. 3n*. 

297.605 

116.77V 

242.26? 

125.10) 

*4092.000 

270.459 

as.  1  TT“ 

7  V.OSo  “ 

2&6.014 

To 0.8 2  V 

253.173 

121. 30-) 

54  128 1 000 "" 

265.106 

23.478 

74.336 

285.644 

121.835 

238.165 

119.2:10 

‘<*155.000 

7r>8,Zrso 

~Tt7F7V 

T 4  TTiU 

2b5.952 

111,962 

246.038 

Tib.  100 

64 1 56.000 

3-7.340 

13.576 

75.672 

284.468 

111.842 

248.158 

1  1  5 .  1  ”  r‘ 

64227.000 

375.587  ' 

2  7  •  7  1 6 

T073T6 

289.684 

1 0  3 . 8  5  2 

256. 148 

1 14.  no 

54263.000 

322.  132 

24.647 

74.032 

285.968 

105.918 

254.082 

113.700 

34294.000 

‘'111:???'" 

24.247 

74  767/ 

285.391 

10. • 54? . 

267.458 

~n  1 .000 

54330.000 

335.491 

25.099 

76,617 

284.383 

102.959 

267.041 

110.80* 

343/1 .000 

541.913 

25.43J 

7o  •  6>  / 

2 &3# 395 

l6l  .492 

2  58.508 

10  4.  no 

543*6. rton 

345.361 

26.770 

78. 174 

281.826 

102.744 

267.256 

1  Co. 90  0 

54429.000 

. 35T.177 

.75771  t- 

- BT57775 

279.241 

lei.rzir 

258.460 

I  06. 000 

54463.50?! 

334.739 

27.536 

80.284 

279.716 

100.777 

259.223 

105.900 

54496.000 

*  356.905 

29.466 

6  71  *  7 1 

292.250 

8  •  26^ 

271*731 

1  06'.  77b? 

"54T32.000 

354.986 

27.094 

74.106 

285.894 

91.807 

268.193 

105.600 

54563.000 

354,679 

4?  1  •  H 

r.  J .  4  5  9 

279.541 

‘  72.56/ 

267.433 

>  <  ...  •  i 

349.967 

17.173 

87.237 

272.763 

97.572 

262.428 

102.700 

34637. 000 

49.085 

3  1  & • 4^ i 

T576T5 — 
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106.200 

107.700 

5573Y52T 

5624.527 

Tinn 

.177 

33.406 

31.099 

91.305 

87.490 

26b.  {>9* 

272.510 

116*637 

112*057 

”7577357“ 

247.943 

"  166.9oo 
105.700 

76Tt 

5687 

T9F7 

.987 

3*5.613 

J7TT85 

28.032 

71^20? 

“7577575“ 

288.793 

106.919 

98.123 

293*0^ 1 
261.877 

lo7.9o* 

105.600 

57275.7T94 

5751.494 

358.905 

359.365 

26.398 

87  75TB 
96.846 

“77V7?57“ 

263.154 

i Q2  *620 
119*118 

“75773 8E“ 
240.882 

■“103.360 

103.300 

57  83790  " 
5614.913 

'  C;iTF?! — 

6.4  24 

2  V  .  5  7V 

28.118 

116.727 

119.261 

“7777771 

240.739 

139.35B' 

138.662 

220.642 

221.338 

162.706 

100.800 

5510 

5941 

•  951 

T.474 

rrnmt 

17TT57T — 

"7357577“ 

145,844 

214.156 

160.366 

”5974 

.351 

r 


»• 


6,401 

30.212 

1  9  7 . 0  H- 

236.982 

144.735 

215.265 

101.900 

6005.351 

“a. *>79 

729.605 

84. *0! 

275.709 

63.601 

296.399 

70.300 

6031.434 

3  1  7,  3  "" 

29.482 

3?.  'it. 

•>27,264 

9  V.  2  1 6 

264.784 

1  2 e- .  6  "  0 

60-62. 4  34 

9.179 

26.562 

ir9.(-4-‘ 

231.157 

144,410 

215.590 

98.700 

6  10."  ,  640 

7 . 3t  5 

2  /  .  0 1 0 

— TT-.6M'. 

231.316 

IhS.SW 

214. 044 

99.30' 

0 1  31 .84" 

6.249 

26.334 

125.42k 

231.572 

145.313 

214.687 

100.000 

6163.959 

77947 

mtoi 

D".o26 

■w-ma 

MKHiHM 

mwrmm 

(.194.9  5  > 

5.517 

22.511 

130.21 : 

229. 7e7 

143.258 

216.742 

98.400 

6227.276 

Trrrr~ 

21.989 

131.283 

228.717 

216.469 

98.000 

6£5b*27£ 

6.738 

20.912 

133  .123 

228.877 

141.184 

218.816 

97.000 

6  2  9  r»  •  7  0  4 

HHnKBIKM 

19.417 

129.207“ 

230.193 

rTTToce 

222.994 

95.500 

63.1. 704 

2.493  • 

18.564 

126. e9a 

233.402 

134.102 

225.898 

97.800 

6417.802 

6.615 

17.478 

125.00? 

224,498 

13". 024 

2  ?9 ,976 

95.600 

644  S • 6  0? 

5.815 

16.511 

124.047 

236.053 

128.176 

231. R24 

96.20" 

6481.456 

1  ;.7S7 

17.727 

■rams 

KlftlM 

237.518 

93.900 

6512.456 

:  7.655 

17.723 

1HWU 

241.182 

91 .0C0 

•6545.0 1 7 

Is.lTT- 

1  b  .)£  3 

1W.0/4 

242.925 

117.145 

242.656 

92 ,  9(Jo 

6i!  7b. 0  1  7 

18.C06 

18.117 

114.414 

245.686 

114.468 

245.532 

91 .500 

6608.420 

18.267 

18.541 

112.245 

247.755 

112.280 

247.720 

90.700 

6634.429 

18. 11C 

18.269 

106.142 

253.858 

106.220 

253.780 

92.200 

6671.682 

1 7.817 

17.873 

101. 832 

258. 168 

101.860 

258.140 

90.600 

6702.682 

17.898 

17.951 

101.470 

258.530 

101.495 

258.505 

0  C  •  6  C 

6735.263 

17.120 

1 7.049 

94.412 

265.588 

94.376 

.  265.624 

88.700 

6766.263 

343.335 

343.322 

01.069 

268.931 

263.925 

69.4C0 

6758.633 

343.127 

343.196 

90. 77C 

269.630 

90.336 

269.664 

90. 6C  0 

6525.633 

341.142 

340.086 

"8.694 

261.406 

98.661 

261.339 

87.900 

6861.979 

19.826 

14.268 

94.651 

265.449 

92.960 

'  267.040 

87.200 

6392.979 

19.120 

14.758 

*)  *■  •  0  9  ? 

264,917 

94.050 

265.950 

87.»0.' 

691)5.404 

20.476 

12. '736 

96. 7  89 

264,.iil 

. ‘74,351 

£63.619 

&  6 , 1 0  0 

6956 • 404 

19.533 

13.97* 

92.61  7 

91.828 

268.172 

87.200 

69bd. 762 

19.6C3 

i  •  4  l  c? 

94.404 

266.656 

93.715 

266 .265 

6  6  *  4  *j  Cj 

7019.762 

209.626 

2S6.699 

6.802 

353.198 

8.795 

351.205 

88.900 

7054.255 

Ta'i7T92 

250.944 

6.918 

353.082 

6.430 

353.570 

90. 30" 

7065.255 

15.633 

13.192 

94. '101 

266,099 

95.019 

264. 981 

88.800 

:115.49c 

16.308 

10.633 

92.83? 

2  6  7.1*67 

93.461 

266.549 

67.000 

7146.498 

18.235 

11.283 

83.70| 

ntnxMdi 

84.488 

275.512 

86.600 

7175.127 

18.633 

7.198 

34. "8 7 

275.913 

86,611 

273.489 

84.60r 

7210.127 

16.224 

8.515 

83.473 

276.527. 

85.461 

274.639 

86.400 

7242.829 

~Te.646 

5.706 

P0.766 

279,  i1}  5 

84,649 

278.451 

84.00" 

7273.829 

17.270 

9.789 

75.311 

284,489 

78.458 

281.542 

85.90" 

13.941 

—mm 

Vt  •  ^  J>  71 

(6.653 

283. 147 

64.500 

7337.576 

18. *68 

5.984 

71.227 

288. 773 

76.476 

283.524 

84.100 

7370.309 

1 8.157 

10.254 

StVioS 

'  293. boo" 

69.556 

290.444 

86. 5oo 

74oT. 309 

24.207 

5.335 

56.176 

301.824 

67.229 

292.771 

81,400 

7434.103 

757777 

77777 

571737? 

TT57T7B 

61.8.19 

76'S.  101 

7465 . 103 

25.769 

9.924 

47.279 

312.721 

55.187 

304.813 

82.hOC 

7497.938 

11.735 — 

4  2  »■->■*  i 

■<  w.nr? . 

50.902 

■aitrUlOiM 

BTTTr- 

■  II  III  1  111 

26.075 

14.282 

37.753 

322.247 

44,238 

315.762 

84.100 

7561.761 

Jr. 542 

— T577FT 

24. 9  0 4 

339,4*6 

33.651 

326. 349 

3  2'.  Tr.o 

T572 .7FT 

32.215 

16.0*0 

11.76  5 

344,245 

27.787 

332.213 

81.  ICO 

7625., 497 

-  327777 

i 7 • H  84 

*14#  !!  4 

>45,976 

26.182 

■"TTTrrn? — 

BT.vTf- 

7 1*  6  C  •  4  0  7 

29.803 

17.884 

14.954 

345.146 

25.369 

334.631 

63.300 

7666.930 

7772771 

i V • 8  0& 

J  3.64V 

'46. Jjv 

<!b.3s5 

3  14  •  tjl  5 

82.400 

7  *  1  9  .  v  3  0 

32.4C3 

18.015 

12.19? 

'347.808 

26.507 

333.493 

61.400 

(752.193 

J  3  •  V  ?  C> 

1 V • 085 

'  11'.  646 

34b. 3ih 

— 7TT5T9 — 

3  J l • 5b 1 

TSVTCTT 

/ 7  fi  3 . 1 "3 

35.572 

16.260 

12.224 

347.776 

30.911 

329.089 

77.600 

7515.474 

■■■■EfcZXLUB 

■I  CTVNM 

■  H1WU 

'4.4.  fV3 

— 377T78 — 

3  /  .  o  n  2  “ 

/  1.  ,'Or 

7  R4cj  #4  74 

30.975 

363.763 

mmsm 

272.673 

111.981 

248.019 

72.  ICO 

5  5  9 1  .632 

<!  9  •  t’  i  3 

3  34  •  i  l  £ 

t  ’H  •  /'  i  i 

o  9  •  1 1 1 

240 . 040 

1  u 

24.839 

354.127 

94,170 

265.630 

115.940 

244,060 

75.i"- 

9654.455 

7777779 — 

““T5T7TV? — 

7  0  •  n/ 

't  — 

1 i  i  «  n  t4 

l  4  it  •  ''  ?  t * 

(  6  •  r>  .  - 

h  '•  •  v  6 

b 


•9  1 


12. 60S  35"*. 689 


.309  3b'. 786 


114.743 


12S.B30 


2 

4 


324.681 
029.550 
'J  32 « 280 
■>36.071 


.252  23.635 


258.969 


3 

359.984 


35T.?fl7 
355. e24 
•35T.T76 
355 .587 

‘SiTJVT 

353.134 


34.150 


37.446 

TTTTTTo' 
. 36.753 
39.208 
38.847 


121.000 


1 

115.960 


1  1 

JJ n  •JLL'L 

~To'9  .Tn*5 
109.1S0 
loV.2vr" 
109. C35 


38.920 


229. 149 
244,040 


249,001 

250.195 

250.820 

252.710 

250.965 


350.738 


351.325 

352.414 

•35G.7T69 

354.210 

-354-.T124- 

353.537 


43.791 


44.929 

45.074 

44.715 

44.290 

TTTFTT 

44.416 


.50 
115.006 


113.16 
125. 777 
TTZTbir 

126.456 

T?r.nnr 

125.422 


244.994 


241,840 

234.223 

T7T7T7T 

233,544 


234.578 


120.464  239.536 

122.031  237.169 
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120.464 

122.031 

239.536 
237.  169 

77.700 
78. 100 

IHI  ■ 

"  234.673 

231.663 

757400 

81.800 

9767. 166 

9618.166 

228.514 

224.330 

64.700 

86.500 

9650.665 

96  r- 1.665 

13b. 447 
36.283 

221.553 

323.717 

88.200 

89.400 

9914.632 

9945.632 

36.204 

34.261 

323.796 

325,739 

90.900 

93.800 

9975.266 
l  ~  r  9  •  2  (t  fc 

32.489 

29.966 

327.511 
-  330.034 

9*  .90-? 
96  •  3  r>c 

1 oc 41.. 561 

ICC  72.o61 

137,701 

133.438 

,  222.219 
226.562 

97.600 

99.600 

10105.413 

1  0  1  *6 . 4  1 3 

127.423 

124,099 

232,577 

235.101 

102.500 
101. 2 CO 

10169.11? 

10200.113 

“  1777617 
121.715 

2 3  i  •  306 

238.286 

162".‘bi30 

1C4.000 

K-2?2.$5e 

10263.928 

115.591 

117.728 

241.409 

242.272 

1 C  5 . 4  0  0 
107.700 

10296.777 

10327.777 

lib. 347 
117.547 

241.653 

242.453 

106.300 

107.200 

10  360.833 

1 .0  391.633 

117.614 

117.690 

242.386 

242.310 

109. ICO 
11C. 600 

1  '■'624  •  668 
10455.668 

119.662 

120.885 

240.318 

239.115 

110.700 

109.900 

]  “  i  o  8 «  5  2 1 
10519.521 

121.079 

122.918 

*  238.921 
237.082 

112.600 

112.200 

1.0  5  52.547 
10503. 547 

126,585 

127.282 

233.415 

232.718 

lU.000 

115.800 

limit™ 

152. oeT 

134.104 

?2  7. 9 19 

225.896 

114.200 

114.000 

1C66 0,718 
10711.718 

137.061 

138.739 

7277539“ 

221.261 

115.400 

116.100 

107*»4.665 

10775.865 

"  143,560 
149,439 

216.440 

210.561 

116.700 

117.000 

10609.045 

10640.045 

152. lb? 
153,895 

20 1  •  S 1 7 
206,105 

115.300 

115.800 

T0673. 140 
10904. 140 

154.849 
'  156.254 

ToTrm 

203.746 

115. 700 

1 16.200 

T0937TT?  3 
10968.273 

. 

■ 

• 

• 

t 

>, 

l 

1 

T 

i* 

, 

("•I  #  « 


SID  0 ?32 
STCP 
SALL 
SIRJOB 

SIPFTC  EGHEVA 
CEGHEVA 
C 

C  PROGRAM  SHUVE.  NEFDS  FCTRAL.  FIELD 

C  FPOCH  CORRECTED  FIELD  VS.  ALT. 

DIMENSION  A( 1000) «G( 30.301 .Hi 10 .30 ) .FMT ( 1? ) 

1. TGI  30. 10) .THIS0.30) «GCORP(30.3O) .HCORR 130.^0) 

COMMON  HCORB  •  GCORR 
C  KTIMF  .3,  GET  POSITION 

C  TIME  a  TIME  DIFFERENCES  IN  YEARS 

C  KTIME»1 .2  CODFFS  ARE  FOR  MAIN  FIFLO.  PATE  OF  CHANGE 

LOUNT«50  ,  ' 

PI«3. 14159265 
PIDEG«5729.57795E-2 

1  READ  IS. 9000JKPR06. NDUMMY. COMULT, KTH  .KTIMF.TIME 

9000  FORMAT (2I4.E16.8.214.E16.R) 

KOMNTaO 

R2»0.0 

GO  TO  I12.12.2010).KTIMF 
C 

C  NDUMMY  is  NO.  OF  COEFFIf IFNTS. 

C  COMULT  IS  A  MULTUPLIFR  FOR  COEFFS. 

C  tCPROG  IS  1.  A  IS  NOT  MODIFIED.  KPROG  =2.  A«A*COMULT.  ! 

C  KPROG  a  3  FOR  SCHMIDT  A.  KPROG  s  4  FOR  VESTINF  TYPE  IN*SCHMIOT) 

C 

12  DO  112  K>1.1000 

112  A  IK  )>0.0 

C 

C  COMPUTE  N1.M22.M1,  FOR  FBASICi  NTOP  IS  NO.  OF  TERMS  IN  COMPLETE 

C  SET  OF  HIGHEST  DEGREE.  NEXTRA  IS  NO.  OF  TERMS  IN  INCOMPLETE  SET. 

C 

OOOOFL-SORT I FLGA  V  ( ND‘JMMY+1 )  ♦  •  01 ) 

npart-oooofl+i.o 

C 

C  READ  IN  COEFFS.  A10, AH  , Pll  ,A20, *.•••.  . 

C 

9  READ  15.100) IFMTII) »l«l»l?) 

100  FORMAT  1 12A6 ) 

READ  (5. FMT) (AIJ) ,J«1, NDUMMY) 

C  FIND  MODIFIED  COEFFS. 

C 

GO  T0I17. 92.91 .91). KPROG 

92  D093  LL*1. NDUMMY * 

93  A(LL)aA(LL)*COMULT 
GO  TO  17 

91  !«0 

N«0 

AA«2.0*(SQRT I  FLOAT I NDUMMY) )+l«0)+1.0 
KK«AA 

D010000  K>3,KK,2 
NaN+1 

D010000  J-l.K 
I-I+l 

M»II-  N**2  +l)/2 
-IF  I W I  1002*15.16 

15  FACTOR  aFCTRAL  I 2»N) / ( 2.0»«N*IFCTRAL IN)) *»2 ) 

GO  TO  14 


MARCOTTE  EGHFVA  F3 

TIMF=02.PAGFS=1S 

CONTINUE 

OLOGIC 

LIST.RFF.PFCK.SDD 
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1 6  FACTOR  *FCT9AL(2*N>/(2.0**N*FCTRAL(N) >*.50RT(?.n/(FCTRAL(N*M)*  FCTR 
1ALIN-M)) ) 

14  GO  TO  (1007. 1010. 1400**14002) .XPROG 

14003  FACTOR  ^FACTOR  •  COMULT 
GO  TO  14001 

14002  FACTOR  *  FACTOR  *  COMULT  /FLOAT (N) 

14001  A  (  IIsFACTOR  *AMl 

1FM-NDUMMV)  10000.17.1  111 
10000  CONTINUE 
17  IX  sO 

no  2000  NX*? .NPAPT 
DO  20P0  MX*1«NX 
IX-IX+1 

IF  I MX-1) 2001 .2002.2003 
200?  GO  TO  (50.51  I .XT IMF 

50  G(NX.1)*A( IX)  4 

GO  TO  7000 

51  TG(NX .1  )  =  A( IX) 

GO  TO  7000 

7003  GO  TO  (52.531. XTIMF 
53  Tfi ( NX .MX ) * A ( IX) 

IX*IX+1 

TH(NX.MX)«A( IX) 

GO  TO  7000 

52  G ( NX .MX ) =A( I X ) 

IX«IX+1 

H(NX.MX)*A( IX)  I 

2000  CONTINUE  | 

GO  TO  1 

2010  DO  55  JJ*1.30 
DO  55  KK*1 .30 

HCORR(JJ.XX)«H(JJ,XX)+TIMF*tH(JJ.XX) 

55  GCOPR ( JJ.KX ) «G( JJ «X< ) 4T IME*TGI JJ.KK ) 

12010  READ  (5.7000)  CIME.AUNCH.TMETA. PHI .HEIGHT 
7000  FORMAT  ( 2F10.4.3 ( 2XF8.4 ) > 

HEIGHT«HEIGHT*1. 85325 
IFICIME-99999. 0)6001 .4007.7001 
4002  COUNT  *K0UNT 

SSGRxSQRT (R2/COUNT) 

WRITE  (6.5005)SS0R 
5005  FORMAT (E16.C ) 

GO  TO  22 

6001  CALL  FIELDJTHETA. PHI. HEIGHT. NPART.X.Y.Z.F) 

FOUNT  *K0UNT+1 
HFLL»SORT(X#X+Y*Y) 

ANC«PIDEG»ATAN(Ae5(Z/HELL) »  • 

IF(X)500. 501.502  * 

500  D«S!GN(PI.Y)-SIGN|ATAN(Y/X».Y> 

GO  TO  503 

501  D«SIGN(PI**5»Y) 

GO  TO  503 

502  D«ATAN(Y/X) 

503  D«PIDEG»D 

IF (LOUNT-50)  1012.1011.1012 

1011  WRITE(6.1013) 

LOUNT«0 

1012  WPITF(f.llO)  THETA.PHI .HEIGHT. X.V.Z.HFLl  .P.D.ANC.AUNCH 
PUNCH  111 .THFTA.PHI .HEIGHT .X.Y.Z.F .AUNCW 

110  FORMAT(2FR.4.1PE16.«.“?P«F17.3.0P7F17.3.F17.4) 

111  FORMAT! 3F10.3.-2P4F10. 3 .0PF10.3 ) 

LOUNTeLOUNT+1 
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GO  TO 

12010 

100? 

CALL 

DUMP 

1003 

CALL 

DUMP 

1007 

CALL 

DUMP 

1010 

CALL 

DUMP 

ini 

CALL 

DUMP 

2001 

CALL  DUMP 

1013  FORMAT I 1 3?H1  F.  LAT.  N  LONG.  HFIGHT  KM.  NORTH  X  FAST  V 

1  DOWN  Z  HOP  I Z  TOTAL  INTENS.  DECLINATION  INCLINATION  T  A 

2FTEP  LAUNCH) 

22  CALL  EXIT 
STOP 
END 

SIPFTC  IFELO  LIST. REF. DEC*. SOD  <* 

FURROW I NE  F 1  ELD  I OL AT .DLONG .HGT .NMAX »RN .AF.RV.R) 

C  FARTHS  MAGNETIC  FIELD  USING  ANY  SFT  OF  COEFFICIENTS 

DIMENSION  HI  30.30) »G 1 30. 3D) .PI  30.30 1 *DP 130.3ft) . CONST  1 30. 30 ) .SP 1 30) 
1 .CP  1 30) *AOR 1 30 ) 

COMMON  H.G 
!FICPI1)-1*0H.2.1 

1  Pii.n-i.o 

DP  1 1 • 1 ) «0«0 

SPI1)«0.0 

CPI1)«1.0 

00  A  M.1.30  I 

DO  3N-1.2  S 

3  CONST IN»M)»0*0 

00  A  N.3.30 
FM«M 
FN«N 

A  C0NSTIN.M)«|  |FN-2*0)»IFN-7*0)-(.FM-1.0)» IFM-1.0))/I lFN*FN)-3«0) / 1 1 F 

1N+FN ) -5*0) 

2  PHI-DLONG/57.2957795 
AR»6371. 2/16371. 2+HGT) 

C-SIN  I OL AT/57 *2957795 ) 

S-SORT  I1.0-C*CJ 

SPI?) «SIN  |PH!) 

CPI2)«C0S  I PM I ) 

AOR 1 1 ) >  AR*AR 

AOR I 2 ) «  AR*AOR ( 1 ) 

00  5  M«3 .NMAX 

SP|M)«SPI2)»CP|M-1)+CP|2)»SP|M-1) 

CP I M ) «CP 1 2 ) #CP I M- 1 ) -SP 1 2 ) *5P I M- 1 ) 

5  AORIM).  AR*A0R|M-1) 

RV«0. 

RN«0.0  * 

RPHI«0.0 
00  6  N«2 .NMAX 
FN«N 

SUMR«0.0 
SUMT«0.0 
SUMP >0*0 
00  7  M«1 *N 
IFIN-MlB.9.8 

9  PIN.N) «S*P|N-1 *N-1 ) 

DPIN.N)«S»0PIN-1.N-1 )+c»pin-i,n-ii 
GO  TO  10 

8  P IN. M) «C»P IN-1 «M) -CONST IN.M)»P|N-2»M) 

0«»  IN.M)  «C«OP  IN-1  .»»)-S*P  IN-1. M)-CONST  IN.  M)#0P|N-2.M) 

10  FM«M-) 

TS«GIN.M)*CP|M)4H|N.M)*SP|M) 
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- $UMR«$UMR+P»N.M|*TS 
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Figure  27  77 

LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  A 
OUTPUT  VOLTAGE  1 
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LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  A 
OUTPUT  VOLTAGE  2 
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LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  B 
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Figure  30  no 

LEAST  SQUAP.ES  APPROXIMATION  FOR  SUN  SENSOR  B 
OUTPUT  VOLTAGE  2 
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Figure  Si  si 

LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  C 
OUTPUT  VOLTAGE  1 
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Figure  32 

LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  C 
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Figure  33 


LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  D 
OUTPUT  VOLTACE  1 


0s-4S. 795+3. 296V+ 12 . 528V2- 1 .970V3 


40*- 


20#- 


V) 

uj 

a 

a 

o 


-20°- 


t 


! 


1 

4 


V 


VOLTS 


DEGREES 


PI  0  1 


Figure  34 
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LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  D 
OUTPUT  VOLTAGE  2 
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Figure  35 

LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  E 
OUTPUT  VOLTACE  1 
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Figure  36 

LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  E 
OUTPUT  VOLTAGE  2 

S -446. 387-7. 424V-9.191V2*!. 400V3 
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Figure  37 


LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  F 
OUTPUT  VOLTAGE  1 
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Figure  38 
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LEAST  SQUARES  APPROXIMATION  FOR  SUN  SENSOR  F 
OUTPUT  VOLTAGE  2 

0--46.021+3.228V+12. 277V2-!. 895V3 
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APPENDIX  F 

In  this  section  we  will  prove  the  equivalence  of  relations  (66)  and  (69). 
It  has  been  shown  that 


Aj  cosp  -  sinp  *  cosys 


If  we  now  make  the  substitution  by  letting 


tan*  *  j- 
A1 


A  ^  ♦  A  ^ 
A1  A2 


V  2  ♦  A  2 
A1  A2 


Substituting  (73)  into  (72)  yields 


A2“  (cos*  cosp  -  sin*  sinp)  -  cosys 


/  A,2  +" 


J  ♦  Aj  ,  COS  (p  ♦  *)  ■  COSys 


cosy 

cos  (p  ♦  *)  -  rajg- 


C05V 

->  P  »  •  arcos  (  -  * 


:°5Y 

y  .  i  arco5  (5IH-£  ) 


sin  B.  -  cos  y. 

s  s 


sin  Bs  -  cos  ys 


y- 


j 


1 


I 


It  follows  that 


cosy 

'  •'  eos(*v)  -  rr-5- 


>/ sin26s  -  cos2'* 


sin(-v) 


’  ;  Ah? 


8,  -  COS  ,s 


Let  us  restrict  ourselves  to 


*  arcos  C  ) 


(It  can  be  shown  that  if,we  take  the  positive  sign  for  v,  the  analysis  will 
result  in  an  erroneous  expression  for  (66)).  , 


Therefore  from  (74), 


i  i 

/ 


AjCOSYs  Ve^S)  *' 

*  — — -  '  -  ■  i  ■■■ 


sin  6.  sin  B. 

3  3 


Multiplying  and  dividing  by  cos0. 


A1  C0Svs  GOs0  +  c(*#ers) 
- - - r 

sin  B_  cose 
s 


(75)  : 


(76) 


/However 


i  i 


A.  cosy. 1 C°sO  w  cosY«[sin0'  -  sin0,cos8«] 

X  )  5  S3  3 


,so  the  expression  in  (76)  is  in  agreement' with  (66). 


i  i 


i  i 


r- 


i 


102 


In  a  similar  manner  we  find 


/ - 2  2 

♦  A.  /sin  e  -  cos  Y  -  A,  cosY 

5in<) - i - _Ti - » 


or 


sine.  sin*  8. 


-  A.  £***PS)  -  A2  cosy 

sinp  -  - - - »-ry- . — - 5.  (77) 

sin  8. 


The  problem  now  is  to  show  that  sinp  in  (77)  divided  by  cosp  in  (75)  will  yield 
the  same  expression  as  (69). 

•in.  -*l<Vr«  '  *2  “*». 

Multiplying  and  dividing  by  the  conjugate  of  the  radical  .in  the  denominator  we  get 


tanp  •[;  A ^  ^in  6#-cos  ys  -A2cosys]  [AjCOSYj*  A^  ^in^fi'-cos2ys] 
[AjCOSYj  ;  A2/sin^8s  -coaZTs)  [AjCOSYj  1  A2  ^  sir.2 0s '  cos ^ y3 ] 


taip  »  - 


2  .  ! 


A1A2cos2yj  ?  A22cos>s  *4in  g#-cos  Yj  ;Aj2cosYj  ^in^j-cos^-AjAjtein^-cos*^  )| 

I 


■v *"S - 5 - 5 - 5 — 

COS  y j (Aj  ♦  Aj  )  -  Aj  sin^es 


tanp  »  ♦  cosy#  / si"*8,  *  cos  Ys  (A)2  ♦  A22)  -  AjA2  sin20s 


_ — — - y 

sinz6s(cos*Ys  -  A2z) 


tano  •  -AjAj  ♦  cosys  /  sin28s  -  cos  ys 

COS^Y j  -  A22 

or 


- 2 - 2 

s  s*n  ®s  "  cos  Vs 
A22  -  COS2  Yj 

which  is  in  agreement  with  (69). 


( 


